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ABSTRACT 

In this second paper of the series, we pursue the analysis of the 180 ks XMM-Newton 
campaign towards the young open cluster NGC 6231 and we focus on its rich OB 
star population. We present a literature-based census of the OB stars in the field of 
view with more than one hundred objects, among which 30% can be associated with 
an X-ray source. All the O-type stars are detected in the X-ray domain as soft and 
reasonably strong emitters. In the 0.5-10.0 keV band, their X-ray luminosities scale 
with their bolometric luminosities as log Lx — log Lboi = — 6.912±0.153. Such a scaling 

■ law holds in the soft (0.5-1.0 keV) and intermediate (1.0-2.5 keV) bands but breaks 

down in the hard band. While the two colliding wind binaries in our sample clearly 

\Q , deviate from this scheme, the remaining O-type objects show a very limited dispersion 

■ (40% or 20% according to whether 'cool' dwarfs are included or not), much smaller 

than that obtained from previous studies. At our detection threshold and with our 
sample, the sole identified mechanism that produces significant modulations in the O 
star X-ray emission is related to wind interaction. We thus propose that the intrinsic 
X-ray emission of non-peculiar O-type stars can be considered as constant for a given 
star. In addition, the level of X-ray emission is accurately related to the star luminosity 
or, equivalently, to its wind properties. 

Among B-type stars, the detection rate is only about 25% in the sub-type range B0-B4 
and remains mostly uniform throughout the different sub-populations while it drops 
k> significantly at later sub-types. The associated X-ray spectra are harder than those 

of O-type stars. Our analysis points towards the detected emission being associated 
' with a physical (in a multiple system) PMS companion. However, we still observe 

a correlation between the bolometric luminosity of the B stars and the measured 
X-ray luminosity. The best fit power law in the 0.5-10.0 keV band yields logLx = 
0.22(±0.06) logL bol + 22.8(±2.4). 

Key words: Stars: fundamental parameters - Stars: early-type - X-rays: individuals: 
NGC 6231 - X-rays: stars - Open clusters and associations: individual: NGC 6231 
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1 INTRODUCTION 

X-ray emission from early-type stars of spectral type O was, 
in December 197 8, one of the earliest findings of the EIN - 
STEIN satellite (lHarnden et al]|l979l : ISeward et alJll979t) . 
It was soon realized that all O-type stars were X-ray emit- 
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ters. Most of them were characterized as soft (kT < 1 keV) 
and rea sonably strong (10 31 < L x < 10 33 ergs" 1 ) 
sources. lHarnden et alJ (Il979l) already suggested that the 
X-ray luminosity was directly linked to the characteristic 
luminosity of the emitter, though at the time the au- 
thors proposed a scaling law between Lx and the visual 
luminosity. However as the stars of their sample had 
mostly the same colors, it is equivalent to assert that Lx 
scales with Lboi- Using various sam ples of stars observed 
with EINSTEIN, different authors (lLong fe White! Il980|: 
I Pallavicini et al.l Il98ll ICassinelli et all Il98ll : IVaiana et alJ 
Il98ll ISeward fc Chlebowskil Il98^) confirmed the so- 
called canonical relation Lx ~ 10" 7 Lboi- Based on the 
EINSTEIN X-ray Observatory catalog of O-type stars 
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jChlebowski et aUll989h . ISciortino et all dl99(tl performed 
a more comprehensive study of the relationship between 
the optical and X-ray properties of O-type stars. They 
confirmed the existence of a canonical Lx — Lbol relation, 
though with a rather large dispersion. They were unable 
to find any significant correlation with the rotation rate 
(usini), the wind terminal velocity (Vao) or the mass-loss 
rate (M), but observed a strong correlation with the 
wind momentum (Mvoo) and with the w i nd lu minosity 
(0.5Mv% o ). More recently, IBerghofer et alJ Jl997l) investi- 
gated the properties of the bri ght OB-type stars det ected 
in the ROSAT all-sky survey feerghofer et alJll99rjl and 
found that the canonical relation extends down to spectral 
type B1-B1.5. They established the separation line be- 
tween the O star relation and a less constrained relation for 
B stars to lie at a bolometric luminosity Lboi ~ 10 38 ergs -1 . 

Historically, two physically different models have been 
proposed to explain the hot O-type star X-ray emission. 
While no standard dynamo-driven surface magnetic field 
is theoretically ex p ected for any star hotter than A7, 
ICassinelli fc Olsonl l)l979T) suggested a scaled-up version 
of a solar-type coronal emission model that yielded a 
roughly correct prediction for the X-ray flux but could 
not explain the softness of the observed spectrum. Indeed 
in the coronal emission model, the X-rays are produced 
near the photosphere and are thus expected to suffer 
absorption by the overlying dense wind. Such an absorption 
is however not seen in the X-ray spectrum of hot stars. 
This suggests that the X-ray emission is rather produced 
throughout a significant fraction of the wind volume, which 
has been a strong argument in favor of the embedded 
wind-shock model. In this second model, the X-ray emission 
is supposed to arise from shocks occurring in the denser 
layers of the winds and that grow from small-scale instabil- 
ities of t he line - driven winds. Since the phenomenological 
model o f iLucvl Jl982T) . hydrodynamical s imulations (e.g. 
lOwocki et alJ 119881 iFeldmeier et alJ Il997h have brought 
further support to this interpretation. 

However, recent high resolution spectral observations 
of early-type stars are bringing the wind shock model to 
its limits. The unprecedented spectral resolution reached 
by the XMM-iVeuiton and Chandra observatories (corre- 
sponding to about 300 kms -1 in velocity space) now allows 
us to probe the widths and profiles of the X-ray emission 
lines seen in the spectra of hot stars. These contribute to 
put constraints on the velocity (Doppler b roadening) and 
location (fir line ratio, IPorauet et alJl200 ll) of the emitting 
plasma, yielding thus an unprecedented c haracterization 
of its localization in the expanding winds. iKramer et alJ 
(2003) reported that the O-type supergiant £ Puppis 
(04 Ief) displays broad, blue-shifted and asymmetric line 
profiles that are generally consistent with the wind-shock 
model. However, observations of other early-type stars 
suggest different pictures and hybrid magnetic wind models 
have been proposed. For example, the X-ray emission lines 
in the spectrum of r Scorpii (BO. 2 V) are significantly 
narrow er than expected from the standard wind-shock 
model dCohen et aljEoo!? ). These auth ors rather suggested 
magn etically confined wind shocks ljud-Doula fc Owockil 
2002), eventually coupled with the clump infall model 



jHowk et al]l200Cft . as the prime origin for the observed 
X-ray emission. A similar model was also su ccessfully ap- 
plied to the ma gnetic rotator 9 1 Ori C (Q5.5V. |Gagne et alJ 
120051) . Finally. TSchulz et al.1 <2003h have suggested that 
young massive stars could enter the main sequence carrying 
a significant residual magnetic field. Therefore, massive 
ZAMS stars could generate their X-ray luminosities via the 
standard model and magnetic confinement may provide an 
additional source of X-rays. 

Compared to single stars of the same spectral type, 
close massive bin aries are also known to displa y an extra 
X-ray emission llChlebowski fc Garmanvl ^991) which is 
generally attributed to a wind-wind collision. Within the 
interaction region, the shocked gas is expected to be heated 
to temperatures of a few 10 7 K and to generate a substantial 
amount of X-rays, which are thus produced in addition 
to the intrinsic emission by each of the components. This 
extra-emission can be further modulated because e.g. of a 
variation of the optical depth along the line of sight towards 
the wind interaction zone due to the orbital motion. It 
could also reflect the changing properties of the shocks due, 
for example, to a variation of the distance between the two 
stars in an eccentric binary system. 

While the earliest works based on EINSTEIN data 
suggested the canonical relation to ext end through the 
B-type range and down to A5 stars (e.g. iPallavicini et al] 
ll98ll) . ISchmitt et alJ <ll985f) showed that the latter law does 
not hold for A-type s tars. Later works llRosner et al.lll98^ : 
ICassinelli et al.l Il994h further suggested that the X-ray 
emission from B stars of spectral type B2 or later was not 
following the same scheme as O-type stars. Because these 
B-type stars do not have the convective zones required to 
sustain a magnetic dynamo, coronal emission is not ex- 
pected. Their stellar winds are also much weaker than those 
of their hotter sisters, and are therefore not supposed to 
provid e a sizable amount of X-rays. Indeed. IBerghofer et alJ 
( 1997) reported that, for stars of spectral type B2 or later, 
the detection rate drops below 10%. Actually the intrinsic 
emission from B-type stars could, comparatively, be much 
lower than from O -type stars, with log (Lx/iboi) ~ —8.5 
JCohen et alJll997l) . This suggests that most of the detected 
X-ray emission for B stars in distant clusters is actually 
associated with unresolved companions, either in a binary 
system or located by coincidence on the same line of sight. 
The question of the intrinsic B-type X-ray emission has 
however not yet received a satisfactory answer. One of 
the difficulties is the intrinsic lack of homogeneity of the 
B-type population that contains different kinds of objects 
(/3Cep, Be stars, shell stars, ...). An additional difficulty 
is their lower emission level (if any), which thus limits the 
number of detections on shorter duration exposures or in 
distant fields. It was hoped that the advent of the 'large' 
X-ray observatories, which combine improved sensitivity 
and spatial resolution, would help to solve this questio n. A 
scan of the recent literature (e.g. IStelzer et al.ll2003l and 
references therein) indeed favours the 'companion' scenario 
but provides by no mean a definitive answer to this question. 

The young op en cluster NGC 6231 (age ~ 3 to 5 Myr, 
iBaume et aljll9991) is considered as the core of the Sco OB 1 
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association and contains a large OB star population. Lo- 
cated a t about 1.6 kpc ( DM — 11.07 ± 0.04, see discus- 
sion in ISana et alJ l2006d . hereafter Paper I), it offers an 
excellent opportunity to probe a homogeneous sample of 
early-type stars in terms of e.g. distance, age, reddening, 
environment and chemical compo sition. O ur XMM- Newton 
campaign has been described in iPaper j and has revealed 
hundreds of point-like sources in the 15' radius field of view 
(FOV) of the satellite. In the present paper, we focus on the 
O- and B-type star population in the FOV. The analysis 
of the sources associated with optically faint counterparts 
is postponed to a forthcoming paper in this series (Sana et 
al. - Paper III, in prep aration) . Pre l iminar y results of this 
work were presented in ISana et alJ (2006a) but should be 
considered as supplanted by the present analysis. 

This paper is organised as follows. The next section pro- 
vides details of the data handling. Sect. |3] makes a census of 
the OB star population in the core of NGC 6231 and iden- 
tifies the early- type X-ray emitters. Sect, ^investigates the 
X-ray properties of the detected O and B-type stars and 
the appropriate Lx — £boi relations are derived in Sect. 14.21 
Sect. |5] presents the properties of the individual early-type 
sources while Sect. HJ discusses the results of our study. Fi- 
nally, Sect. [7| summaries our main results. 



2 X-RAY DATA HANDLING 

The present work relies on the X-ray source catalogue and 
on the reduced X-ray data as obtained and described in Pa- 
per I. Based on a census of the OB population in the FOV, 
we have now identified the early-type X-ray emitters in the 
FOV (see Sect. |HJ. Using the SAS v6.0 software, we have 
then extracted the OB star X-ray spectra from the merged 
event lists, i.e. from the event lists built using the combi- 
nation of the six pointings. For each object, we adopted a 
circular extraction region centered on the X-ray source po- 
sitions listed in Paper I and with a radius corresponding to 
half the distance to the nearest neighbouring X-ray source. 
Due to the crowded nature of the cluster core in the X-rays 
(see e.g. Fig. 1) and to the limited spatial resolution of the 
EPIC detectors, the background could not be evaluated in 
the immediate vicinity of the stars, but had to be taken from 
the very few source free regions. We adopted several circu- 
lar background regions spread throughout the FOV. For the 
two EPIC MOS instruments, as the background can be con- 
sidered as uniform on a particular CCD, we associated, to 
each X-ray emitter, source-free background regions located 
on the same CCD detector. Regarding the pn instrument, 
the background is rather considered to depend on the dis- 
tance to the detector read-out nodes. We therefore selected 
source-free regions approximately situated at the same dis- 
tance from the read-out nodes as the target of interest. The 
adopted background extraction regions for MOS and pn in- 
struments are displayed in Fig. Their respective positions 
and sizes are given in Table Q We used the appropriate re- 
distribution matrix files provided by the Science Operations 
Centre (SOC), according to the position of the considered 
source on the detectors. We built the corresponding ancil- 
lary response files using the arfgen command of the SAS 
software. The spectra were finally binned to have at least 10 



Figure 1. From top to bottom, combined EPIC MOS1, MOS2 
and pn images in the range 0.5-10.0 keV. The background extrac- 
tion regions adopted according to the different positions of the 
sources on the detectors have been overplotted. North is up, East 
to the left. A colour version of these figures is available in the 
electronic version of the paper. 
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Table 1. Positions and sizes (r e xtr.) of the adopted background 
extraction regions in the FOV of the different EPIC instruments. 



Bkg. 


Label 


a (J2000.0) 


8 (J2000.0) 


f*cxtr. 


MOS1 


CCD 1 


la 


16 h 54 m 31!43 


_41°45'42'.'2 


20" 




lb 


16 h 54 m 23?28 


-41°46'14'.'7 


20" 




lc 


1 fi h 53 m 44 s 1 2 


—41 ° 53'34"fi 

rt-L <JiJ J1. \J 


25" 


CCD 2 


2a 


16 h 55 m 12!05 


-41°49'34'/0 


50" 




2 b 


16 h 54 m 57?95 


-41°49'16'/7 


50" 




2c 


16 h 54 m 52?31 


-41°44'13'.'0 


35" 




2d 


16 h 54 m 58!19 


-41°41'29'.'2 


45" 


CCD 3 


3a 


16 h 54 m 38?42 


-41°39'33'/3 


30" 




3b 


16 h 54 m 18?46 


-41°39'36'.'0 


40" 




3c 


16 h 54 m ll!88 


-41°36'51'.'0 


30" 


CCD 4 


4a 


16 h 53 m 28?04 


-41°40'06'.'8 


35" 




4b 


16 h 53 m 13?61 


-41°42'10'.'2 


45" 




4c 


16 h 53 m 22!95 


-41°45'35('4 


45" 


CCD 5 


5a 


16 h 53 m 04?39 


-41°49'59'.'9 


50" 




5b 


16 h 53 m 21!91 


-41°56'55'.'4 


60" 


CCD 6 


6a 


16 h 53 m 47?44 


-42°01'OtX'9 


45" 




6b 


16 h 54 m 12!79 


-42°02'23'/5 


45" 


CCD 7 


7a 


16 h 54 m 56?95 


-41°58'55'.'4 


50" 




7b 


16 h 55 m 09f87 


-41°53'17'.'6 


60" 


MOS2 


CCD 1 


la 


16 h 54 m 31!43 


-41°45'42'.'2 


20" 




lb 


16 h 54 m 23?28 


-41°46'14'/7 


20" 




lc 


16 h 53 m 44!12 


-41°53'34'.'6 


25" 


CCD 2 


2a 


16 h 54 m 12!79 


-42°02'38'/5 


45" 




2 b 


16 h 54 m 55!16 


-41°58'50'/5 


50" 


CCD 3 


3a 


16 h 55 m 09?87 


-41°53'17'.'6 


60" 




3b 


16 h 55 m 12!05 


-41°49'35'.'0 


50" 




3c 


16 h 54 m 57?95 


_41°49'16'.'7 


50" 


CCD 4 


4a 


16 h 54 m 52!53 


-41°43'55'.'5 


30" 




4b 


16 h 54 m 57?29 


-41°41'14'.'2 


60" 




4c 


16 h 54 m 38?53 


-41°39'29'.'6 


30" 




4d 


16 h 54 m 24?48 


-41°39'38'.'5 


25" 


CCD 5 


5 a 


16 h 53 m 49!69 


-41°37'34'.'7 


40" 




5b 


16 h 53 m 28?04 


-41°40 / 06'/8 


35" 


CCD 6 


6a 


16 h 53 m 22?95 


-41°45'35'.'4 


45" 




6b 


16 h 53 m 01?20 


-41°49'59'.'8 


50" 


CCD 7 


7a 


16 h 53 m 21!91 


-41°56'55'.'4 


60" 




7b 


16 h 53 m 47?44 


-42°0lW9 


45" 



counts per bin. The analysis of these spectra is presented in 
Sect-H 

Finally, holding the source positions fixed at the val- 
ues given in iPaperl we ran the SAS task emldetect on the 
individual pointings to obtain the background-subtracted, 
vignetting- and exposure-corrected count rates of the ob- 
jects during each observation. Similarly, we used the extrac- 
tion regions described above to build X-ray good time in- 
terval corrected and background subtra cted ligh t curves in 
the different energy bands considered in lPaperl The latter 
curves were computed for a series of time-bin sizes ranging 
from a few hundred seconds to 5 ks. These temporal infor- 
mation allow us to probe the variability of the early-type 
X-ray emitters on the different time-scales covered by the 
present data set (see Sect.Hj. 



Table 1 - continued 



Bkg. 


Label 


a (J2000.0) 


<5 (J2000.0) 


r e xtr. 


pn 


A 


8 b 


16 h 53 m 51?25 


— 41°37'04'.'7 


30 




7b 


16 h 54 m llf 43 


— 41°36'48'.'5 


30" 


B 


9b 


16 h 53 m 22f 96 


— 41°39'56'/7 


30 




7a 


16 h 54 m 20f 24 


— 41°39'36'.'0 


30 




10b 


16 h 54 m 38?42 


— 41°39'33'.'3 


30 




lie 


16 h 54 m 58?19 


— 41°41'29{'2 


45" 


C 


10a 


16 h 54 m 31?43 


— 41°45'42"2 


on// 

20 




lib 


16 h 54 m 54?54 


— 41°44'28yo 


35" 




12b 


16 h 55 m 24?85 


— 41°46'43"3 


60 


D 


9a 


16 h 53 m 22?95 


-41°45'35'/4 


45" 




8a 


16 h 53 m 42!61 


-41°45'38'.'3 


30" 




11a 


16 h 55 m 04?20 


_41°48'31'.'5 


35" 




12a 


16 h 55 m 14!27 


_41°48'57'.'4 


30" 


E 


6a 


16 h 53 m 12?40 


-41°49'02'.'6 


30" 




6b 


16 h 53 m 23?00 


-41°50'44'.'2 


30" 




5 a 


16 h 53 m 40?35 


-41°49'14'.'5 


25" 




2a 


16 h 55 m 01?01 


-41°52'06'.'6 


35" 




3a 


16 h 55 m 14?79 


-41°52'49'.'9 


60" 


F 


6d 


16 h 53 m 23?29 


-41°53'2E'2 


40" 




5b 


16 h 53 m 39?86 


-41°55'02'.'0 


20" 


G 


6c 


16 h 53 m 20?12 


-41°56'32'/9 


40" 




4a 


16 h 53 m 54?85 


-41°58'36'.'0 


30" 




la 


16 h 54 m 28?93 


-42°00'10 / /9 


45" 




2b 


16 h 54 m 54?93 


-41°58'53'.'0 


50" 




3b 


16 h 55 m ll?85 


-41°58'18('7 


45" 



3 EARLY-TYPE STARS IN NGC 6231 

We first take a census of the OB star population in 
the observed FOV. Relyin g on a large review of the 
literature (see refere nces i n iP aper fl) . on the Catalog of 
Galactic OB Stars jReedl 1200311 . and on the WEBDA 1 
and SIMBAD 2 databases, we selected all stars for which 
at least one of the previous references quoted a spectral 
type O or B, regardless of their sub-spectral type, their 
luminosity classification or their single/multiple status. 
This resulted in 106 objects among which 91 B stars and 
15 O stars. These are listed in Tab le 2 av a ilable from the 
CDS (http://cdsweb.u-strasbg.fr). iHoukl <197Sh quoted 
HD 152437 as a B9 st ar. However , we recently revised its 
classification to A0III JSanall2005l) and estimated its dis- 
tance to about 500 pc. This object is thus not included in the 
present census. Together with various cross-identifications 
(Cols. 1 to 7), Table 2 provides the adopted positions (Col. 
8) and the V and B — V magnitudes (Cols. 9 and 10) 
obtained from the S SB06 cat alogue (Sung et al. 2006, in 
preparation; see also iPaper J for a brief description) . The 
visual (Mv) and bolometric (Mbol) absolute magnitudes 
(Cols. 11 and 12) were computed using a distance modulus 
DM = 11.07, the intrinsic c o lours and bolometric correc- 
tion scale of ISchmidt-Kalerl il982T) and a red dening law 
A v = R x E(B - V) with R = 3.3 iSung et alJll99gT) . A v 
was computed for each object to acc ount for the d ifferential 
reddening existing across the cluster jSung et alll998Tl . The 
adopted spectral types and the corresponding reference are 

1 http://obswww.unige.ch/webda/ 

2 http:/ /simbad. u-strasbg.fr/Simbad/ 
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Table 3. List of the early-type X-ray emitters in the XMM-JVeiuion FOV. The first five columns give the source identifiers. Col. 6 
indicates the separation d cc between the associated optical and X-ray sources. Cols. 7 and 8 respectively provide the adopted instrument 
combination and the radius of the source extraction region. In Col. 7, the notations MI and M2 respectively stand for MOSI and MOS2. 
The adopted background regions are indicated in Cols. 9 to 11 for each EPIC instrument (see Table0. The last three columns provide 
the bolometric luminosity (Lbol)i the ISM-corrected X-ray luminosity (Lx) an <i the Lx/^bol ratio respectively. In the present table, the 
objects have been ordered by decreasing Lbol- 







Object IDs 








Instr. 


f*cxtr . 


Backj 


;round rej 


;ions 


log L hol 


log L x 


L 

log 7i 




HD 


CPD-41° 


Se68 


SBL98 


(") 




(") 


Ml 


M2 


pn 


(ergs- 1 ) 


(ergs" 1 ) 




[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


279 


1 52248 


7728 


291 


856 


0.8 


EPIC 


19.0 


I 


I 


t; 


39.49 


32.87 


—6.62 


202 


1 52234 


7716 


290 


855 


1.2 


MOS 


10.0 


I 


I 




39.45 


32.42 


—7.03 


216 


1 ^99'^ 


7718 


306 


858 


1 3 


EPIC 


14.0 


I 


1 


]j 


39.30 


32.40 


—6.90 


292 


152249 


7731 


293 


857 


1.4 


EPIC 


11.5 


I 


I 




39.22 


32.30 


—6.92 


291 


152247 


7732 


321 




5 


EPIC 


23.0 




5 


jj 


38.99 


32.09 


—6.90 


422 


326331 


7744 


338 


571 


0.7 


MOS 


12.0 


I 


I 




38.92 


31.96 


—6.96 


468 


152314 


7749 


161 


615 


8 


EPIC 


15.5 


I 


1 


]j 


38.90 


31.80 


—7.10 


185 


1 ^991 8 


7713 


2 


853 


1.1 


MOS 


20.0 


3 


5 




38.77 


31.87 


—6.90 


149 


1 ^991 Q 


7707 


254 


234 


9 


EPIC 


08.0 


I 


I 




38.68 


31.80 


—6.88 


306 




7733 


297 


862 


1.1 


EPIC 


11.5 


I 


1 




38.62 


31.83 


—6.79 


372 




7742 


224 


505 


1.9 


MOS 


13.0 


I 


I 




38.48 


31.94 


—6.54 


126 


1 52200 


7702 


266 


206 


0.7 


EPIC 


10.0 


I 


I 




38.44 


31.29 


—7.15 


6 


152076 


7684 






1.2 


MOS 


37.5 


4 


6 




38.44 


31.04 


-7.40 


566 a 




7760AB 


810 




1.6 


pn 


17.0 






G 


38.44 


30.19 


-8.25 


313 


326329 
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1.0 


EPIC 


10.0 
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36.08 


30.73 
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265 


200 


1.2 


EPIC 


09.0 


1 


1 


E 


35.64 


30.81 


-4.83 



a See individual comments in Sect. 15.21 



quoted in Cols. 13 and 14. For binaries, the quoted visual 
and bolometric magnitudes are given for the system as a 
whole. We used the respective visual brightness ratio of the 
system components and derived a reddening correction that 
accounts for the component spectral types and for their 
relative contribution to the total light from the system. 
Doing so, we deduced the My for the system. We then 
derived the individual visual and bolometric magnitudes. 
Finally, we combined the latter values to obtain the total 
bolometric magnitude for the system. The last column 
(Col. 15) of Table 2 provides miscellaneous notes about the 
objects, together with the relate d references . Comp ared to 
the previous results presented in ISana et alJ (I2006afl . some 
spectral types have been revised thanks to recent FEROS 
observations <ISanal l2005). The distribution of the selected 
stars among the different spectral su b-types and luminosity 
classes given in ISana et alJ <l2006al) still provides a good 
overview of the early-type population in NGC6231. 

As a next step, we cross-correlated the X-ray catalogue 



of |Paper"3 with the obtained list of OB stars in the FOV. We 
first built the distribution of the number of associated OB 
counterparts as a function of the cross-correlation radius. 
The curve reveals 30 associations within a cross-correlation 
distance (d cc ) of 2'.'5. Six additional sources are found with 
dec between 3'.'0 and 3'.'6. The remaining association s require 
dec > 4'.' 8. For consistency with the results of IPaper fl. 
we adopted a cut-off radius of 2'.'5, thus resulting in 30 
OB X-ray emitters (see Table 01. Except in one case, the 
associated early-type stars are the only optical counterpart 
found in the SSB06 catalogue within the adopted cut-off 
radius. Considering that faint optical sources located close 
to an early-type star (which is intrinsically bright) are 
difficult to detect and might indeed be missing in the 
SSB06 catalogue, one can however estimate the probability 
to spuriously associate an OB star with the X-ray emission 
produced by a lower mass star loca ted by chance on the 
same line of sight. Using the lPaper J results and considering 
only stars with V < 19, one obtains that this probability 
remains below 0.03. As a consequence, the early-type stars 
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Table 4. Best-fit single temperature models of the form: wabsisM * wabs * mekal. The first and second columns give the X-ray and 
optical identifiers. Col. 3 provides the adopted interstellar column (N^ M , in 10 22 cm -2 ). The next three columns (Cols. 4 to 6) list the 

best-fit absorbing column (JVjj, in 10 22 cm -2 ), temperature (kT, in keV) and normalisation factor (in cm -5 , norm = ^ d a 

■ / N c N H dV 

with d, the distance to the source — in cm — , N B and JVjj, the electron and hydrogen number densities - in cm -3 ). The quoted upper 
and lower values indicate the limits of the 90% confidence interval (see text). The reduced chi-square and the number of degrees of 
freedom (dof) of the fit are given in Col. 7. Cols. 8 and 9 provide the observed (fx) an d the ISM absorption-corrected (/£ orr ') fluxes (in 
10 -14 ergcm — 2 s -1 ) in the range 0.5-10.0keV. The last column provides the sum of the numbers of source counts in the fitted spectra. 
The horizontal line separates the O-type objects (upper part of the table) from the B-type stars (bottom part of the table). 



[1] 


ID 

[2] 


N XSM 
[3] 


7V H 
[4] 


kT 

[5] 


norm 
[6] 


xl (dof) 
[7] 


fx 
[8] 


f corr . 

•fx 

[9] 


Net cts 
[10] 


6 


152076 


0.29 


< 0.12 


0.28°;* 


1.85f Jg x 10 -5 


1.08 (92) 


0.92 


3.42 


1095 


22 


326343 


0.36 


43 - 87 


o-568:l| 


< 4.73 x 10 -5 


1.08 (13) 


0.25 


0.57 


170 


167 
448 


326320 
7746 


0.24 
0.30 


61 73 


61 ' 82 

U - D1 0.50 

n 58 - 75 

U.OSq 4 3 


2.50 3 ;|| x 10 -5 

2 - 99 B:88 x 10-5 


0.94 (70) 
1.12 (11) 


0.79 
1.01 


1.24 
1.92 


762 
144 


566 


7760AB 


0.42 


7Q 1 - 69 


°- 28 ao8 


< 1.71 


0.76 (23) 


0.18 


0.48 


265 



listed in Table [3] are most probably physically related to 
the observed X-ray sources. 



28 are detected, which constitutes a much smaller detection 
rate. 



After a visual inspection of the X-ray images at the 
position of the undetected OB stars, we further added two 
objects to the list of early-type X-ray emitters. Indeed it 
clearly appears that the related locations on the EPIC 
images have been corrupted because of the crowdedness 
of the field in these particular regions. The spectra of 
these sources could of course not be securely extracted 
because of the probable contamination by very close X-ray 
neighbours. These two sources are thus not included in the 
discussions of Sect. 0] still they are very probable X-ray 
emitters. These two additional sources are identified with a 
cross (+) appended to the ID number of the X-ray source 
(Col. 6) in Table 2. For these two objects, we emphasize 
that th e corresponding X-ray source in the catalogue of 
IPaper J is probably constituted of different components 
that the psf (point spread function) fit has not been able to 
separate. Finally, we found residuals larger than expected 
for source-free regions at the locations of seven B-type 
stars. These might be X-ray sources be low the detection 
threshold of our catalogue (see details in lPaper m . For the 
sake of completeness, we also identified these stars in Table 
2. 



As a first result, the 15 O-type objects in the FOV 
are all seen in the X-rays while only ~25% (15 out of 63) 
of the early B stars (sp. type B0-B4) are detected. This 
percentage remains almost constant within the subtypes 
B0 to B4 but reaches ~ 45% among the giants and su- 
pergiants (4 out of 9). Except for the respectively known 
(CPD-41°7737) and suspected (Se 209) slowly pulsating 
B-type stars (SPB) in the FOV, that are both detected in 
the X-ray domain, no particular correlation is found for the 
different sub-populations. Three out of the six known /3Cep 
are detected (alternatively, three out of the eight known or 
suspected /3 Cep are detected). Five out of the 16 known 
binaries are also detected (alternatively, six out of the 26 
known or suspected binaries are also detected). For the later 
B-type stars (subtypes B5 and later), only two objects out of 



4 X-RAY PROPERTIES 
4.1 X-ray spectrum 

Using the procedure described in Sect. |5J we extracted the 
merged X-ray spectra of each OB emitters. We selected the 
deepest instrument combination (see Table^ , requiring that 
the source is not located too close to a gap or a detector edge. 
The radius of the source extraction region and the identifi- 
cation of the associated background regions in the different 
instruments are given in Table [3] The spectra are relatively 
soft and peak at an energy well below 1.5 keV. These X- 
ray spectra are presented in Figs. 2 to 5, also available from 
the CDS. Beyond the different emission levels, the O- and 
B-type spectra show clearly different characteristics. 

To quantify the physical properties of the X-ray emis- 
sion, we adjust ed a series of optically thin therm al plasma 
mekal models ijMewe et alJIlQSflt lKaastralll992T l to the ob- 
tained spectra. For the purpose of scientific analysis, we lim- 
ited the considered energy range to the 0.5-10.0 keV band 
and we simultaneously fitted the di fferent EPIC s pectra us- 
ing the XSPEC software v.11.2.0 3 l)Arnaudl|l996i ). Because 
of the expected presence of local absorption (e.g. due to 
the winds), we requested, in the spectral fits, a total equiv- 
alent column density of neutral hydrogen larger or equal 
to the interstellar (ISM) column density N^ M . The latter 
was estimated from the color excess using the relation of 
iBohlin et ail l)l97Sh : 7V H = E(B - V) x 5.8 x 10 21 cm -2 . 
Fits with a common local absorption column for multi- 
temperature models were, for various cases (e.g. HD 152234, 
HD 152248, HD 152249, ...), unable to reproduce the X- 
ray spectra. We thus adopted independent local absorb- 
ing columns for each mekal component in the model. As 
no strong deviation from the solar metallicity is seen in 
the optical spectra, we kept the chemical composition fixed 
at the solar values. We successively tested one-, two- and 
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Table 5. Same as Table|^for two-temperature models of the form: wabsjgM * (wabsi * mekali + wabS2 * iekal2 ) . normi and norm2 
are expressed in 10 cm - J . 
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/ IUI 1 1 b\ 
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[2] 




[3] 


W 


[5] 
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[7] 


[8] 
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oo24.70 
- oi 7.37 


402'22 
"'^"0.30 


7lS Zf 

"■ ' ^0.65 


7 992 22 
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.22 


222' 2 2 
"•"0.12 


27S- 2 . 9 

"■^'0.25 


105 


"•"'0.58 


"•'"0.72 


35 62 9 § 
J "-"31.2 


1.47 (234) 


39.74 


82.23 


12255 


251 


7721p 





.24 


418-1? 

"■^0.31 


249- 28 
"•^0.22 


34 52 4 '2. 

J ^-"30.2 


1 38H2 

±.OOq 4 g 


1 33} -29 


2 48f 29 
^•^°1.44 


0.92 (204) 


4.07 


8.56 


2797 


291 


152247 





.28 


< 0.06 


25S-2? 

"•^"0.21 


1 4922-2 
^•^13.0 


"■^0.03 


632'92 

"'" J 0.59 


7 4l| f7 

' -^^5.89 


1.26 (351) 


12.20 


37.90 


6932 


292 


152249 





.27 


252 3 ? 

"■^"0.21 


23?- 2 i 

"•^ J 0.22 


145 2 92' 


±,±o 0.86 


8l9'22 

"• Oi 0.73 


12 3i 6 ; 6 


1.40 (389) 


21.39 


61.85 


17380 


372 


7742 





.28 




^qO-62 
"• 09 0.50 


7 48.O 
'• 4 6.2 


7S°- 86 
"• ' J 0.62 


1 05 1 - 26 


10 ql5.9 

ro.» 12 g 


1 HQ /^DQ*! 


14.0 


27.2 


4790 


122 


326331 





.28 


0.29°t| 

0.19 


0.23°-25 

0.20 


77.8" 3 4 

39.7 


0.88 


0.88.V91 
0.76 


lo.e^ 8 

7.9 


1.19 (157) 


10.70 


28.74 


3325 


78 


Se 41 


0. 


.21 


0-65JJ;! 


u - z± 0.19 


2L6 8 1 


76 1,8 
"■ ' D 0.06 


1 66 2 -^6 


L18 a77 


1.11 (74) 


1.31 


2.02 


1198 


121 


Se 265 





.26 


"■ OZ 0.26 


19°-23 


18.7 76 g 7 


n 4Q103 

"■^ a 0.24 


1 83 2 - 47 
i.odj 43 


1 Q7I.88 
1 "-"0.77 


1.37 (105) 


1.34 


2.22 


789 


139 


7706 





.25 


62°- 84 
U - Oz 0.37 


26°- 38 
U - ZD 0.19 


13.2fV 


< 0.51 


cfi3.25 
z -°"1.82 


i.98 2 : 84 


0.99 (92) 


3.02 


4.62 


1183 


203 


7715 





.25 


0.58 ; 2g 


23°- 27 
U - ZJ 0.18 


12.8f 3 3 2 


1 no 2 - 15 


1 32 1 - 82 
1 - oz 1.08 


i.87?: 89 


1.29 (94) 


1.42 


2.32 


1078 


213 


Se 259 





.26 


7^0.87 
U -'°0.61 


n 1 5 016 

U - 10 .14 


185.?!?: 


< 1.77 


2 00 3 - 33 
z,uu 1.26 


1 60 303 

1- DU 1.06 


1.10 (100) 


1.75 


2.88 


1151 


282 


Se 209 





.24 


76°- 83 
U -' D 0.64 


16 017 
U - 1D .14 


169. 280 - 


71 1 ' 05 
(1 0.32 
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three-temperature absorbed mekal models. According to 
the quality of the observed spectrum, we finally adopted the 
model with the lowest number of parameters that still ade- 
quately reproduces the data. Good spectral fits are generally 
obtained using two-temperature (2-T) mekal models with 
independent local absorbing columns for each component. 
However, the brightest sources may require an additional, 
higher temperature component (3-T models) while the spec- 
tra of a few very faint sources are satisfactorily described by 
a single temperature (1-T) component. The observed emis- 
sion is generally soft (kT ^ 2.0 keV) and slightly absorbed 
(A^h ^ 10 22 cm -2 ). Best fit results are given in Tables0]to 
^together with the observed and ISM-absorption corrected 
fluxes in the range 0.5-10.0 keV. Tables [1] to HJ also pro- 
vide the 90% confidence intervals on the best fit parameters. 
These have been computed using the xspec command error. 
In this procedure, each parameter is varied until the value 
of the x 2 ! minimized by allowing all the other non- frozen 
parameters to vary, is equal to the best fit value increased 
by a value Ax 2 = 2.706, equivalent to the 90% confidence 
region for the considered parameter. 

Direct comparison of the spectral properties obtained 
with models composed of a different number of mekal com- 
ponents is often not an easy step. For this reason, the dif- 
ferent result tables were organised according to the number 
of mekal components used in the model. Table probably 
provides the best illustration of the clearly distinct X-ray 
spectral characteristics of the O- and B-type stars. With the 
temperatures of the two mekal components mostly below 
1 keV, the X-ray emission from the O-type stars is partic- 



ularly soft. The first temper ature is about 0.3 keV, except 
for CPD -41° 7742. However. lSana et alJ <2005l) showed that 
the intrinsic X-ray emission from this system is generally 
contaminated by an extra-component resulting from an on- 
going wind interaction process, except when the interaction 
region is occulted. At these particular phases at which we 
probably only detect the intrinsic emission from the O-type 
primary, it is interesting to note that the tem perature of 
the ' cool' component drops to about 0.35 keV flSana et alJ 
2005), thus quite close to the average 'cool' temperature of 
the other O-type objects. 

The absorbing column associated with the lower tem- 
perature component is generally very low. The best-fit value 
might drop to zero so that only an upper limit can be pro- 
vided. It is known, from past experience, that suppressing 
the local absorbing column of the soft component can lead to 
a better fit and to more stable results (see e.g. ISana et alJ 
|2004 l2005h . This suggests that the only effective absorp- 
tion for the related component is due to the ISM. When 
appropriate this is indicated by a 'n.' in Tables HH and H 
The temperature of the second component is usually about 
0.7 keV and the associated column of absorbing matter re- 
mains moderate, with a mean value about 0.5 x 10 cm - . 
We note that the described pattern for the X-ray emission 
from O-type stars as revealed by the 2-T models is also seen 
in the results of the 3-T models (Table The first two 
components are almost unchanged. An additional high tem- 
perature component (with kT > 2 keV) however describes 
the higher energy part of the spectrum. This third compo- 
nent has generally an emission measure at least five times 
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Table 6. Same as Table |4*| for three-temperature models of the form: wabsjgM * ^ (wabSi*mekali ) , with i = 1,3. The fourth column 
gives the number i of the considered component in the model. 
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lower than the less energetic components and is probably 
only seen because the corresponding sources are bright X-ray 
emitters, which allows us to probe the higher energy tail of 
their spectrum. Alternatively, the high energy components 
could come from an unresolved PMS companion located 
along the same line of sight (though, as stated above, this 
is rather unlikely). These objects h ave indeed X-ray spectra 
typic ally harder than O-type stars iFeieelson & Mo ntmerld 
1999), but should not bias the measured X-ray flux of the O 
stars as the latter is at least one order of magnitude brighter. 

The X-ray spectra of the B-type stars show quite differ- 
ent properties (see Tableland Figs. 2 to 5). Compared to 
the O stars, the low energy component has a slightly lower 
temperature but a much larger absorbing column. The main 
difference however is seen for the high temperature compo- 
nent which is systematically larger than 1 keV. Beyond this 
spectral scheme, it is obvious that the B-type X-ray emit- 
ters are fainter than the O-type emitters and, probably for 
this reason, we never needed a 3-T model to satisfyingly 
reproduce their spectrum. 



4.2 Z/x — Lhoi relationship 

To derive the intrinsic X-ray luminosities, we corrected the 
observed fluxes quoted in Tables 2] to |5] for the interstellar 
reddening. We then computed the corresponding unab- 
sorbed fluxes and, adopting DM = 11.07 ± 0.04, the X-ray 
luminosities. Fig.HJpresents the logLx vs. logLboi diagram. 
The different properties of the X-ray emitters (spectral 
type, luminosity class, multiplicity) and their location in the 
log — log plane are given in the upper panel of the figure. 
The lower panel gives, in addition, the upper detection limit 
for the undetected B-type stars, computed according to 
their respective positions on the dete ctor. The se limits were 
estimated using Eq. 5 and Fig. 6 from lPaper J and refer to a 
source with a temperature about 0.7-1.5 keV and a typical 



ISM reddening corresponding to an equivalent column of 
neutral hydrogen of 0.26 x 10 22 cm -2 . We emphasize that 
no additional (local) absorption column is considered in 
this typical model. 

Clearly the O- and B-type stars show different be- 
haviours. The former ones are brighter (as expected) and 
present a steeper Lx — £boi relation than the latter. The 
B-type stars and the late O-type stars also seem to present 
a larger dispersion around the observed trends. However, 
the impression of dispersion towards lower Lx mainly 
comes from five points, namely HD 152076 for the O-type 
stars and CPD-41°7760, CPD-41°7746, HD 326320 
and HD 326343 for the B-type sample. These sources 
correspond to five out of the seven OB sources with an 
off-axis distance over 10'. The other two off-axis sources 
are HD 152235, a Blla star whose X-ray emission could 
not be disentangled from a brighter neighbour (and which 
is thus not plotted in the current diagram; see also Fig. |U] 
in Sect. and HD 152247, an X-ray bright 09 III star. 

With the exception of the latter object, all the studied 
sources with an off-axis angle larger than 10' display a 
smaller X-ray luminosity compared to the general trend in 
the sample. Further inspection of the results of Tables 0] to 
HJ reveals that the X-ray fluxes from these five sources were 
all obtained with single temperature models. These were 
characterized by kT about 0.25 keV for the O star and 
about 0.5-0.6 keV for the four other B stars. Interestingly, 
the obtained kT corresponds to the lower temperature 
component of the 2-T models. We caution that, because 
these sources are located in the outer parts of the detector, 
where the sensitivity significantly decreases, only the softest 
part of the spectrum produces a significant number of 
photons. The harder part of the spectrum provides too 
faint a contribution for a 2-T model being used, which 
probably yields an underestimation of the observed flux in 
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Figure 6. Top panel: ISM-absorption corrected X-ray lumi- 
nosities in the 0.5-10.0 keV band plotted vs. bolometric lumi- 
nosities. The different symbols indicate the different properties 
of the sources. Spectral type: O (filled symbols), B (open sym- 
bols). Luminosity class: supergiant (I or black), giant (111 or red), 
main sequence (V or blue). Multiplicity: binary (triangles), ra- 
dial velocity variable (squares), presumably single star (circles). 
Bottom panel: same as the upper panel, the plain lines give the 
iBerghofer et al.Hl997ft relations obtained in the 0.1-2.0 keV band, 
respectively for objects with logLx > 38 (ergs -1 ) and for those 
below. The dashed line indicates our typical detection limit at a 
distance of about 5' from the axis for a source with a temperature 
kT around 0.7-1.5 keV and an ISM column of absorbing hydrogen 
of AT^ SM = 0.26 X 10 22 cm -2 Csee lPaoer Hh The crosses indicate 
the upper limits for the undetected B-type stars (see text). 
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Figure 7. ISM-absorption corrected X-ray luminosities in the 
0.5-10.0 keV energy band plotted vs. bolometric luminosities. The 
symbols have the same meanings as in Fig. |H] Best-fit linear re- 
lations in the log — log plane for O (Eq. |5J and B (Eq. |7J stars 
are indicated by the dashed and the dotted lines respectively. The 
two vertical bars in the upper part of the graph give the expected 

1- a deviation for B (o"b) and O (co) type stars. 

the 0.5-10.0 keV band. HD 152247 being about one order 
of magnitude brighter in the X-rays, the higher energy 
part of its spectrum remains relatively well constrained. 
The obtained spectrum is indeed of good quality and a 

2- T model is definitively required to properly reproduce 
it. To preserve as much as possible the homogeneity of 
our sample, both in terms of location in the cluster and in 
terms of the hypotheses that underlie the analysis, we have 
restrained our discussion to the results obtained with the 
2-T and 3-T models. 

Before deriving the Lx — Lbol relationships for our sam- 
ple, we would like to draw the attention of the reader to the 
following point. In the rest of this paper, we will make a 
difference between a scaling law (SL): 

Lx = K Lboi, 

where K is the scaling constant, and a power law (PL): 
L x = K Lboi, 

where F is the power law index. In the log — log plane, both 
relations yield a straight line: 

log L x = log Lboi + log K 

for the scaling law, and 

log L x = T log L bo i + log K 

for the power law. We emphasize that the former expression 
has a slope equal to unity in the log — log plane and is equiva- 
lent to the statement log (Lx/Lboi) = est, with est = logK. 
The earliest works on the canonical relation mostly used 
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Figure 8. ISM-absorption corrected X-ray luminosities in different energy bands plotted vs. bolometric luminosities. The symbols have 
the same meanings as in Fig. Best-fit linear relations in the log — log plane are indicated by dotted and dashed lines respectively for 
B- and O-type stars. 



scaling relations while, later on, power laws were used to 
provide an additional degree of freedom to the fitted model. 
In the following, both kinds of models will be tested and 
their results will be compared. 



4-2.1 X-ray emission from O-type stars 
(L bol > 10 38 ergs- 1 ) 

As seen in Figs. |S| and [7] the O-type stars mostly concen- 
trate on a line in the Lx — £boi plane. As a first step in 
the analysis, we computed the linear-correlation coefficient 
r (Table |7J. We then estimated the probability for an un- 
corrected parent distribution to yield a value of r equal o r 
larger than the observed value (see e.g. iBevingtonl Il969f) . 
For all the considered O star samples (see below), we could 
reject this hypothesis at a significance level of 0.001. Assum- 
ing first that Lx/ibol = est, a least-square fit performed on 
the 14 objects (sample A) yields: 

log (L x /iboi) A = -6.865 ± 0.186 (1) 

with a typical dispersion of the X-ray luminosity of about 
55% around Eq. However, we note that most of the dis- 
persion comes from four points. Among them, HD 152248 
and CPD —41° 7742 are known to display a significant extra- 
emission associated with a wind interaction phe nomenon oc- 
curring in both systems l|Sana et, al.l2004l . l2005h . This extra- 
emission, averaged over the orbital cycle may account for at 
least one third of the observed flux (see individual notes 
for more details). Rejecting these two systems from the fit 
(sample B), we now obtain: 

log (Lx/L b oi) B = -6.912 ± 0.153 (2) 

with a typical 1-er deviation corresponding to about 40% 
on L x (see Fig. . Beyond the two colliding wind systems, 
the remaining dispersion is mainly due to HD 152200 and 
HD 326329. It is interesting to note that, with a respective 
spectral type of 09.5 and 09.7, these two presumably sin- 
gle dwarfs are the two coolest O-type stars of our sample. 
Hence, further excluding the stars with a sub-spectral type 
later than 09 (sample C), the best fit values are almost 
unchanged: 



log (L x /Lboi) c = -6.925 ± 0.087 . (3) 

The residuals drop this time to 20%. Henceforth, 
with the exception of the very late O-type main-sequence 
stars and of the two identified systems with strong wind 
interaction effects, the X-ray luminosity of the other objects 
displays a very limited dispersion around the canonical rela- 
tion. Among these ten stars, dwarfs, giants and supergiants 
are found, as well as presumably single stars and binaries 4 . 
In Sect. |S| we pay a more particular attention to each of 
the X-ray emitters. 

As a next step, we computed the X-ray luminosities of 
the sample stars in different energy bands. Plotted against 
the bolometric luminosities (Fig. , we again observed lin- 
ear relations in the log — log plane, both in the soft (Sx: 
0.5-1.0 keV) and medium (M x : 1.0-2.5 keV) bands, while a 
large dispersion is seen at higher energies (H x : 2.5-10.0 keV). 
In the S x and M x bands, the relative position of the stars 
remains very similar. The two colliding wind systems display 
a larger deviation in the intermediate band than in the soft 
band. This is to be expected as the X-ray emission produced 
in the wind-interaction region is typically harder than the in- 
trinsic emission from O-type stars, which drops significantly 
beyond 1-1.5 keV. Using the sample B, best-fit relations in 
the soft (0.5-1.0 keV) and intermediate (1.0-2.5 keV) energy 
ranges yield respectively: 

log (L x ,s/L b oi) B = -7.011 ± 0.193 (4) 
and 

log (L x ,M/iboi) B = -7.578 ± 0.205 (5) 

with a typical dispersion on i x of about 55% and 60% re- 
spectively. The hard band shows a much larger dispersion 
and further suggests a slightly different behaviour for main- 
sequence stars compared to more evolved objects. The for- 

4 The X-ray contribution of a wind interaction to the total emis- 
sion of the latter binaries is expected to be very limited (see 
Sect.|FJ. Their X-ray emission should thus reflect very closely the 
intrinsic emission of their components and, indeed, their X-ray 
luminosities do not deviate from the canonical relation. 
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Table 7. Comparison between best-fit Lx — L bo i scaling law (SL) and power law (PL) models, computed in the 0.5-10.0 keV energy 
bands (upper part of the table) and in the soft (Sx) and intermediate (Mx) energy ranges (bottom part). The first two columns specify 
the O star sample used (see text) and the number N of stars in the sample. The third column provides the log Lx — log L bo i correlation 
coefficient r. The next two columns indicate the best-fit scaling laws, the reduced \ 2 and the number of degrees of freedom (dof) of the 
fit. Cols. 6 and 7 yield the same information for the best-fit power laws. Note that, in Cols. 5 and 7, the x 2 statistics have not been 
normalized by the intrinsic variance of the data. The last column gives the values of the statistics F x . 



Sample 


N 


r 


SL (Eq. #) 


xl (dof) 






PL 




Xl (dof) 




Fx 


A 


14 


0.89 


Eq.IU 


3.46 x 10" 2 


(13) 


log L x 


= 0.91(±0.13)logL boI - 


3.6(±5.1) 


3.62 x 10~ 2 


(12) 


0.4 


B 


12 


0.91 


Eq.H 


2.34 x 10^2 


(11) 


log L x 


= 0.85(±0.12) logL bol - 


1.2(±4.7) 


2.24 x 10~ 2 


(10) 


1.5 


C 


10 


0.97 


Eq.|3 


7.56 x 10~ 3 


(9) 


log L x 


= 0.88(±0.08) logL bol - 


2.2(±3.1) 


6.58 x 10" 3 


(8) 


2.3 


B 


12 


0.86 


Eq.El 


3.75 x 10~ 2 


(11) 


log L x ,s 


= 0.80(±0.15) logL bol + 


0.6(±5.9) 


3.53 x 10~ 2 


(10) 


1.7 


B 


12 


0.85 


Eq.El 


4.22 x 10" 2 


(11) 


log L X ,M 


= 0.75(±0.16)log L bol + 


0.6(±6.2) 


3.95 x 10" 2 


(10) 


1.8 



mer ones are about at least one order of magnitude brighter 
than the giants of a similar spectral type. There also seems 
to be a steep trend in the giant and supergiant population 
and their fluxes in the hard band indeed increase by more 
than one order of magnitude towards larger bolometric lumi- 
nosities. We caution however that the measured X-ray flux 
above 2.5 keV is probably little constrained, except for the 
brightest X-ray emitters, and that the present results for the 
hard band probably require confirmation. 

As mentioned earlier, we also adjusted power laws to 
the data by means of a linear regression performed in the 
log — log plane. The best-fit relations are quoted in Table 
These only yield a marginal improvement of the resid- 
uals and the dispersion of the data around the power law 
relations is virtually unchanged. To test the judiciousness of 
including this additional degree of f reedom in the mod el, we 
used a F^-test as described in e.g. iBevingtonl |l969). The 
F x statistic is given in Table and was obtained from the 
equation: 

2 2 

T? _ Xsl ~ Xpl (( .\ 
X ~xk/(iV-2) (b) 

where N is the number of points in the considered sample. 
Clearly, the F x -test does not indicate any significant 
improvement of quality of the fit, even when adopting a 
significance level as large as 0.1. In consequence, we consider 
that the power law model does not provide any significant 
improvement compared to the simpler scaling law. In the 
following, we have adopted the canonical Lx — L b oi relation 
given by Eq.0 



4-2.2 X-ray emission from B-type stars 
(L bol < 10 38 ergs" 1 ) 

The present sample of B-type stars is strongly dominated by 
non-detections for which only an upper limit on the actual 
X-ray luminosity can be obtained. To account for this, we 
have attempted to app ly various survival analysis techniques 
(e.g. llsobe et aljll986tl . Unfortunately, we could not achieve 
consistent results, probably because of the much larger num- 
ber of upper limits compared to the actual measurements. 
Therefore we have chosen to restrain our analysis to the de- 
tected objects only, making thus the implicit hypothesis that 
the latter stars form a specific population, distinct from the 
undetected B stars. 



As mentioned in Sect. 14.11 the detected B-type stars 
have a quite harder emission (kT2 > 1.0 keV) compared to 
the O-type stars. They are however much more luminous in 
the 0. 5-10.0 keV band than predicted by thelBcrghofc r et alJ 
lll997l) relation (Fig. Ujjl. This could result from the fact that 
the higher energy part of the B-type spectrum falls outside 
the observable energy range of the ROSAT satellite. Like 
the O-type stars, their X-ray and bolometric logarithmic lu- 
minosities are apparently still linked by a linear relation. 
The linear-correlation coefficients r are about 0.75, which 
rejects the null hypothesis of such a high value occurring by 
chance from an uncorrelated distribution at a significance 
level of 0.01. Only the r value computed for the Sx band 
is slightly lower (r = 0.68) and yields the rejection of the 
null hypothesis only if we adopt a significance level of 0.05. 
It is obvious, from Figs. |S]to|5] that the slope of a putative 
relation is significantly different from unity and indeed the 
F x test described above indicates that the power law pro- 
vides a very significant improvement of the quality of the fit 
compared to a simple scaling law. From our selected sample 
of 11 B-type stars, we obtain: 



log L x = 0.22(±0.06) log L bo i + 22.8(±2.4) 



(7) 



with a typical dispersion around Eq.|7|of 0.14 in the log — log 
plane, thus corresponding to about 37% in the X-ray lumi- 
nosities. 

From Fig. the emission of the detected B-type stars 
in the log — log plane also follows a linear relation in the 
different sub-energy bands considered. By opposition to the 
behaviour of O-type stars, this linear behaviour holds in the 
harder band too. However, we note that the softest part 
(< 5 keV) of the hard band is responsible for most of the 
flux. The latter remains indeed very limited above this value 
(see e.g. Fig. 4). From linear regressions, we obtain: 

log L x ,s = 0.21(±0.08) log L bol + 22.6(±2.8) (8) 

log L x ,m = 0.22(±0.07) log L bo i + 22.4(±2.5) (9) 

log Lx,h = 0.24(±0.07) log L bol + 21.1(±2.6) (10) 

with typical dispersions of respectively 0.17 (47%), 0.14 
(39%) and 0.15 (42%) around these relations. Before pre- 
senting a more detailed discussion of the derived Lx — 
L bo i relations and comparing these with earlier works (see 
Sect. the next section will first review the individual 
properties of the early-type X-ray emitters. 
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5 INDIVIDUAL OBJECTS 

This section aims at providing specific information on the 
various early-type X-ray emitters detected in the XMM- 
Newton FOV. In particular, we haverm -R S investigated 
their variability on the different time-scales allowed by our 
data set. Using the emldetect net count rates obtained in 
Sect. |5J we performed a \ 2 test °f hypothesis with a sig- 
nificance level of O.Of and we actually tested the null hy- 
pothesis of a constant count rate throughout the six point- 
ings of the XMM-Newton campaign. Nine sources displayed 
consistent indications for variability in all the EPIC instru- 
ments. These are noted 'var.' after their name in the in- 
dividual sections here below. In addition, CPD —41° 7715 
and CPD -41° 7737 gave consistent results in all the EPIC 
instruments with a significance level of 0.1 only. They are 
noted 'var?'. Finally, HD 326329 presented significant vari- 
ations at the 0.01 level but was probably contaminated by a 
neighbouring strong source (X#324). The latter displayed a 
strong flare that is also seen in the HD 326329 count rates. 
It is thus difficult to determine whether, in addition to the 
contaminating component, the HD 326329 source displays 
intrinsic variation of its flux or not. 

To search for intra-pointing variability (i.e. on time 
scales shorter than the various exposure durations), we fur- 
ther applied the Kolmogor ov-Smirnov and p robability of 
variability tests described in ISana et al.l i2004f) for the case 
of HD 152248. While the latter non-parametric tests do not 
make any a priori hypothesis on the form of the variability, 
they do not allow to correct for the background. We thus 
also built background and GTI corrected light curves us- 
ing various bin sizes ranging from a few tens of seconds to 
a few kiloseconds and we applied similar % 2 tests as previ- 
ously done for HD 152248. As in the case of inter-pointing 
variability, we adopted a significance level of 0.01 for the var- 
ious tests. For the O-type stars, intra-pointing variability is 
only seen for the two strong colliding wind binaries while, for 
B-type stars, the tests reveal several flaring sources. Thanks 
to the amplitude of the flares, the latter sources had how- 
ever already been detected as variable by the inter-pointing 
variability study. 



5.1 O-type stars 

- HD 152076 is located at 13' N-NW of the cluster 
core. It lies outside the pn detector FOV and is thus only 
seen in the two MOS instruments. It presents a constant 
flux throughout the campaign. The quality of the obtained 
spectra is definitively too low to constrain a 2-T model. 
We suspect the adopted single temperature model to bias 
the computed flux in the 0.5-10.0 keV band. This source 
was thus not included in the derivation of the Lx — £boi 
relationship. 

- HD 152200 is the latest O-type star of our sample. It 
is reasonably bright and well isolated in the X-ray images 
and its X-ray flux remains constant. We considered it to 
be reasonably well constrained. HD 152200 is however one 
of the O-type stars that deviate the most from the derived 
Lx — I/boi relation. We have not been able to identify a 
cause that could possibly bias the determination of its 




Figure 9. HD 152235: EPIC three-colour image: red, Sx band; 
green, Mx band; blue, Hx band. The location of X#179, indi- 
cated by the 3" radius circle, is clearly biased by the presence 
of two close, unresolved sources that can only be distinguished 
in this image thanks to their different spectral properties. The 
supergiant B star HD 152235 is clearly associated with the soft 
emission. The image appears in colour in the electronic edition of 
the paper. 

flux and we consider the observed deviation to be significant. 

- HD 152218 displays a slight modulation of its flux 
by about 20% but its X-ray s pectral properties remain 
almost unchanged l|Sanal l2005li . The observed modula- 
tions could be linked to a possible wind-wind collision 
occurring in this O+O system, but the latter authors 
note that second order effects, such as radiative inhibition, 
might actually govern the wind interaction properties in 
this system. However, the HD 152218 X-ray flux remains 
probably dominated by the intrinsic emission from the stars. 

- HD 152219 (var.) presents a slight modulation of its 
flux jSana et alJl2006bfl which could be related to a similar 
wind interaction as in CPD —41° 7742. Except for the longer 
orbital period of HD 152219, the two systems are indeed 
very similar. The amplitude of the variation seems however 
much more limited than in CPD— 41° 7742. This probably 
results from the longer period of the system, yielding thus a 
larger separation between the components and a subsequent 
dilution of the wind material that enters the interaction. 
HD 152219 is the third O-typ e eclipsing binary known in 
the cluster JSana et alJl2006Hl 

- HD 152233 is most probably a long period binary (P ~ 
months. ISanal 12005 V but its orbital properties are still not 
constrained. The HD 152233 X-ray flux presents an appar- 
ently constant level throughout the campaign. 

- HD 152234 (var.) is reminiscent of HD 152233. The 
system is mo st probably a long period (P ~ months) binary 
(Sana 2005), but its orbital properties are poorly con- 
strained so far. Its X-ray flux slightly decreased throughout 
the XMM-iVeuiton campaign with an average slope of 
about —4.3 x 10 -3 cuts -1 per day. This corresponds to a 
reduction from about 75 x 10 -3 to about 55 x 10 -3 cnts -1 
between the first and the sixth observations. This could be 
due to a wind interaction, between the two components, 
whose strength is varying because e.g. of the changing 
separation between the stars. However, a more precise 
analysis requires first to constrain the orbital properties 
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of the system. As for HD 152233, an ongoing monitoring 
campaign should soon provide this missing piece of the 
puzzle. 

- HD 152235 (not resolved) is a Blla star in the south- 
ern part of the field. The cross-correlation did not associate 
this object with an X-ray source. The three-colour image 
(Fig. |5J however clearly reveals that HD 152235 is a very 
soft X-ray emitter. The source detection has actually been 
biased because of a close, harder emitter. Due to this close 
neighbour, the HD 152235 spectrum can unfortunately not 
be reliably extracted and the source was thus not included 
in the present analysis. 

- HD 152247 is located in the outer regio n of the field. It 
is most probably an SB1 binary l)Sanal l2005). but its period 
is largely undetermined. It could however be of the order 
of one year. The X-ray flux seems to remain constant from 
one pointing to the other as well as within each pointing. 

- HD 152248 (var.): the colliding wind binary HD 152248 
is the brightest X-ray source in the FOV l|Sana e t al. 2004) 
and presents phase-locked modulations of its X-ray flux due 
to the varying strength of the shock in this eccen tric system. 
The X -ray data have been thoroughly analysed in lSana et alJ 
(2004). To provide averaged properties for this peculiar sys- 
tem, we adjusted a 3-T model to its spectrum. Unfortu- 
nately, we could not get xt values better than 2.4. This 
might reflect the changing spectral properties of the sys- 
tem and the impossibility for mekal models to reproduce, 
in this case, an averaged emission level. Because of the x2 
> 2, the XSPEC software did not allow to explore the pa- 
rameter space to search for a better local minimum. For 
the same reason, we could not directly compute the limits 
of the 90% confidence intervals. In Table HJ the quoted in- 
tervals were estimated from the 1-u error bars that were 
converted into 90% confidence intervals assuming a Gaus- 
sian distribution of the uncertainties. We note however that 
the obtained flux in the 0.5-10.0 keV band is very similar to 
the avera ge of the fluxes ob served during the six individual 
pointings iSana et alJl2004T) . This indicates that the present 
result is quite consistent with the previous analysis. Using 
the newly derived Lx — £boi relation of Eq.|5|and the indi- 
vidual luminosities, we expect log Lx = 32.58. We observed, 
on av erage, log Lx = 32. 87, which is about twice larger. 

In 1 ana et al] J2004D . we have modelled the X-ray 
emission expected from the wind-wind interaction using 
hydrodynamical simulations. Once diluted by the intrinsic 
emission from the two stars, the modelling qualitatively 
reproduced the observed modulations. However, the pre- 
dicted X-ray flux was about a factor two too faint compared 
to the observed emission level of the system. At the time, 
we attribute d this difference to th e observed dispersion 
around the iBerghofer et al] <ll997tV s canonical relation. 
From the present results, this dispersion is presumably 
much more limited than previously expected. Therefore, 
the HD 152248 system probably deserves further theoretical 
work to quantitatively explain the observed modulations 
and the exact shape of its light curve. 

- HD 152249 presents a strong increase of its flux during 
Obs. 2. This increase is seen in the three energy bands, 



but is most prominent in the intermediate and in the hard 
bands. The corresponding images reveal a nearby flaring 
source located at about 6" S-W of HD 152249 (thus well 
within the extraction region). While a slight contamination 
is still detected in the hard band during Obs. 3, the 
count rates measured during Obs 1, 4, 5 and 6 remain in 
good agreement and no additional intra- or inter-pointing 
variability could be found. Both the HD 152249 spectral fit 
given in Table |S] and the derived fluxes and luminosities 
rely only on the four pointings with no contamination. 

- HD 152314 is a relatively isolated bright X-ray emitter 
whose X-ray flux seems to remain constant. 

- HD 326329 (contaminated?) is located in a relatively 
crowded region of the cluster. The count rates obtained via 
the SAS task emldetect are clearly affected by the bright 
X-ray source X#324 that displayed a flare during Obs. 4. 
We re-derived the count rates in a more limited extraction 
radius of 7'.'6 but obtained mitigated results. Clearly the de- 
rived fluxes for this object should be regarded with caution. 
It is interesting to note that, after the two wind interaction 
systems, HD 326329 presents the largest deviation from 
Eq. |21 From the above remark, this dispersion is probably 
an artefact and might not be linked to any physical process. 

- HD 326331 reveals particularly broad lines dis- 
playing slight p rofile variations in its optical spe c trum . 
Some authors Jffill et all Il974t iLevato fc Morrelll Il983l: 
iGarcfa fc Mermilliodl 1200 ll) have suggested t hat it coul d 
actually be a binary system. From our work dSanal 120051) . 
we believe that the binarity is not confirmed and that the 
profile variations are probably not, or not only, due to the 
presence of a companion. The X-ray flux of HD 326331 
remained constant over the duration of the campaign. 

- CPD — 41° 7721 is a visual double star with com- 
ponents, labelled 'p' and V, separated by ~ 5.8". Other 
denominations for these components are SBL 350 and 351 
or BVF 12 and 27. The X-ray source X#251 is clearly 
associated with the O-type 'p' component. Abnormally 
large residuals are seen at the position of another star, 
SBL 358 (V — 15.89), a neighbouring object, but actually 
there seems to be no detected X-ray emission associated 
with the B-type star CPD -41°7721s. The emldetect count 
rates of the 'p' component show some deviating points, 
especially at Obs. 2 for the MOS2 instrument and at 
Obs. 4 for the pn. These deviating points are however not 
confirmed by the other instruments. We re-derived the 
count rates using the same extraction region as the one 
adopted for the spectra and we found that the X-ray flux 
from CPD — 41°7721p is probably constant throughout the 
campaign. We note that the extracted spectra might be 
contaminated by the residual X-ray emission located at 
~8" from CPD-41°7721p and probably associated with 
SBL 358. The contamination is however probably very low 
and limited to the soft part of the spectrum, thus in the 
energy range where CPD — 41°7721p displays the strongest 
X-ray emission. 

- CPD 41° 7733 presents a n approxim atively constant 
X-ray flux throughout its orbit jSanal 12005'! and indeed, no 
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Figure 10. CPD-41°7706 background corrected EPIC MOS 
light curves (in the range 0.5—10.0 keV) during the XMM.-Newton 
campaign. The average count rates during each pointing are given. 
The vertical bars are the l-tr deviations while the horizontal bars 
give the duration of each pointing. 



significant X-ray overluminosity from a possible wind in- 
teraction is expected to be found at our detection threshold. 

- CPD 41° 7742 (var.) is the second known SB2 
eclipsing earl y-type binary in NGC 6231 and has a period 
close to 2.4 d JSana et alJl2005t . It is abo ut 3 times b r ighter 
in the X-rays than expected from Eq. |5] S ana et alJ (2005) 
presented the MOS X-ray light curve of the system, almost 
fully covering the orbital cycle with a time bin of 1 ks. 
We showed that an extra-emission component is associated 
with the inner face of the secondary. When the orientation 
of the line of sight is favourable, this component could be 
responsible for an increase of the X-ray flux by about a 
factor of 2. These features were interpreted as the primary 
wind crushing on (or nearly on) the secondary surface. 
Alternatively, the wind-photosphere interaction could be 
altered leading to a wind-wind interaction region located 
close to the secondary surface. Even at its lowest flux level, 
the logarithmic X-ray luminosity is still about 0.3 above the 
one expected for the intrinsic emission of this O+B binary. 



5.2 B-type stars 

Before reviewing the individual B-type emitters, we note 
that only about 20% of the B-type sta rs are as sociated 
with an X-ray source in the catalogue of IPaper"! Among 
these, four sources are apparently single and, at our detec- 
tion threshold, do not display significant variations. These 
are Se 127 (= SBL 746, B8III/IV), Se 142 (= SBL 748, 
B4V), Se 160 (= SBL 593, B3III), and HD 326343 (B3III). 
These stars will not be discussed any further. 

- HD 326320 is an SB2 binary located in the outer 
regions of the FOV. At our detection threshold, we do 
not detect any significant variations but the count rates 
are quite low. Its spectrum is well described by a single 
temperature model. 



- CPD-41°7706 (va r.) is a Bl+Bl binar y 
jLevato fc Malarodal Il98d : iGarcfa fc Mermilliodl 1200 if) . 
ISemaTil^OOSl) reported that the Ha and all the Hei lines 
from 14920 to 17100 display mixed absorption plus emis- 
sion. All these lines have the same profile with a narrow 
absorption superposed on a broader emission which results 
in line profiles reminiscent of the Be spectral signature. 
Because the component signatures are blended, it is not 
possible to decide whether both components or just one of 
them is actually a BlVe star. The source presents a smooth 
variation of its X-ray flux (Fig. EOJ, but no sign of flaring 
activities could be found for this system. 

- CPD — 41°771 5 (var? ) is reported as a spectroscopic 
binary by iRab oudj lll996t) and as a f3 Cep-type star by 
lArentoft et alJ feOOlT) . There might be an increase of the 
detected flux during Obs. 2 but a \ 2 test with a significance 
level of 0.01 did not allow us to reject the null hypothesis 
of constant count rate throughout the six pointings. 

- CPD-41°7737 (var?) is a slowly pulsating B2V star 
jArentoft et alJ l200lf) . There is a slight increase (~40%) 
of the count rate during Obs. 5, but because of the large 
error bars, it might not be significant. A \ 2 test with a 
significance level of 0.01 did not allow us to reject the 
null hypothesis of constant count rate throughout the six 
pointings. Note that another optical source, SSB06 #3549 
(V = 16.9), lies at about l'.'9 from X#337, thus within the 
adopted cross-correlation radius. In this paper, we have 
assumed CPD —41° 7737 to be the only physical counterpart 
to X#337. 

- CPD— 41°7746 was r eported as an SB1 system by 
IGarcfa fc Mermilliodl (1200 ll) . Located in the outer part of 
the FOV, it falls on a gap of the pn instrument. The MOS 
count rates are too low to investigate the intra-pointing 
variability. However, no sign of flaring activities is seen. 
The spectrum is very poor and the star was not considered 
further in the present analysis. 

- CPD— 41° 7755 (var.) is an apparently single B1V 
star. The X-ray flux seems in a higher state during the last 
three observations than at the beginning of the campaign. 
An inspection of the images does not reveal any contami- 
nating source. The maximum emission is seen during Obs. 5 
and is about a factor 3 to 4 higher than during the first 
three pointings. We caution however that the count rates 
are low (between 1 and 4 x 10 -4 cnts -1 in the two MOS 
instruments) . 

- CPD 41° 7760 is a visual binary with separation 
of about 2". The A and B components have a probable 
spectral type B0.5III and B1V respectively. Though the 
fainter component (CPD — 41°7760B) lies slightly closer 
to the X-ray source location, the brighter component 
(CPD — 41°7760A) remains within the adopted 275 cross- 
correlation radius. In Table |3] and subsequent figures, the 
adopted L^oi corresponds to the joint contribution of the 
two stars. Lying in the outer region of the field, outside 
the MOS detection area, this source is only seen in the pn 
instrument. Its X-ray emission level is close to our detection 
limit, but seems to remain constant throughout the six 
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Figure 11. Se 41/SBL 149 EPIC pn (top panels), MOS2 (mid- 
dle panels) and MOS1 (lower panels) light curves (in the range 
0.5-10.0 keV) during XMM-Newton Obs. 1 to 3. The time is 
given in ks since the beginning of Obs. 1. The bin size is 1 ks. 
The effective duration of Obs. 2 has been reduced because of a 
high background event. 
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Figure 12. Se 209 / SBL 394 EPIC pn light curve (in the 
range 0.5-10.0 keV) during XMM-Newton Obs. 1 to 6. The time 
is given in ks since the beginning of Obs. 1. The bin size is 2 ks. 
The effective durations of Obs. 2 and 4 have been reduced because 
of a high background event. 

pointings. The obtained spectrum is of poor quality (Fig. 5) 
and was not considered in the derivation of the Lx — L bo i 
relations. 

- Se 41 / SBL 149 (var.) is a B4V star according to 
SIMBAD. With an average pn (vignetting and exposure-time 
corrected) count rate of about 5 x 10~ 3 cnts -1 , it is a 
relatively bright X-ray source. Its count rate is higher by a 



factor 2 to 3 during Obs. 2, whereas it is almost constant 
during the 5 other exposures. This increase is also seen in 
the three energy bands. Inspection of the images indicates 
no sign of a contamination by a nearby flaring source within 
the extraction region. We therefore consider the flare to be 
associated with Se 41. We built light curves with a temporal 
resolution of 1, 2 and 5 ks for all three instruments. They 
consistently indicate that the emission is decreasing from 
the beginning of Obs. 2 and reaches a similar emission level 
as during the other pointings at the end of Obs. 2 (Fig. lllH . 

- Se 209 / SBL 394 (var.) ha s been proposed either as 
a binar y candidate (lRatoud| |l996) or as a slowly pulsating 
B star iArentoft et alJl200ll) . It presents an increase of its 
emission during Obs. 3 (Fig. 1121 , which yields an average 
count rate higher by about a factor 2.5 compared to other 
pointings. 

- Se 259 / SBL 317 is susp ected to display radial 
velocity vari ations llRaboudl Il996l) possibly related to an 
SB2 nature iSanal 120051) . This source seems to present a 
constant, though low, count rate throughout the campaign. 

- Se 265 / SBL 200 (var.) is a late B-type main- 
sequence star (B8.5V accor ding to simbad) and was quoted 
as a cluster non-member bv lBaume et al.l il999T) . Its average 
count rate during Obs. 4 is about a factor two larger than 
during the other pointings while, actually, its emission level 
decreases throughout Obs. 4. However, the count rates are 
quite low (4 x 10 -3 cnts" 1 in the pn and 10 -3 cnts" 1 in 
the MOS, on average during the low emission state) and 
caution should be the rule. 



6 DISCUSSION 

6.1 The O-type X-ray emitters 

In Sect. 4, we derived various Lx — Lboi relations corre- 
sponding to the O-type stars of our sample. The preferred 
relation is the scaling law log(Lx/Lboi) = —6.912 ± 0.153. 
Previous works however proposed power law Lx — Lboi re ~ 
lations and we can use the relations quoted in Table [7] for 
comparison. The latter ones are indeed formally equivalent 
to the scaling relations derived in Eqs. to |3 Working in 
the 0. 2-4.0 keV band with EINSTEIN data. lSciortino et all 
( 1990) proposed 

log L x = 1.08(4-0.06/ - 0.22) logL bol - 9.38(4-2.32/ - 0.83). 

In the 0.1-2.0 keV band. lBerghofer et alJ lll997l) obtained 

logLx = 1.13(±0.10) logLboi - 11.89(±0.38), 

while in the 0.3-12.0 keV band, lAlbacete Colombo et all 
J2003I) reported 

logLx = 1.07(±0.04) logL bol - 6.2(4=0.1) 

(though see Sect. 16.1,211 . Our adopted slope for a power 
law Lx — Lboi relation in the log — log plane is 0.85 ± 0.12 
which is thus s ignific antly lower than the values from 
iBer ghofer et alJ (Il997tl and from lAlbacete Colombo et alJ 
j2TO3Tl ~ Even our scaling relation, with a slope equal to 
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unity, is still at more than 1-cr from the values quoted in 
the two latter works. The identification of the strong CWB 
systems in our sample and their subsequent rejection from 
the fit can explain part of the observed difference. Indeed 
including the CWB systems in the fit yields a larger value 
for the slope (see Table EJ. 

The direct comparison of the present results with 
previous works should however be considered with caution 
because of the different energy bands considered and, in 
particular, of the lower energy cut-off adopted for the 
analysis. The influence of the ISM-absorption correction 
dramatically increases towards lower energy. For example, 
assuming an ISM absorption column of 0.26 x I0 22 cm -2 , 
which is typical for NGC6231, the correction factor is 39.3, 
2.68, 1.43 and 1.04 respectively in the 0.2-0.5, 0.5-1.0, 1.0- 
2.5 and 2.5-10.0 keV bands 5 . Hence, a slight uncertainty on 
the measured fluxes at low energy (or on the estimated ISM 
absorption) will have a large influence on the intrinsic fluxes 
computed. The very lo w energy edge (0.1 keV) adopted by 
IBerghofer et alJ il997T) could be one of the reasons for the 
large dispersion observed in their sample. Another reason 
is pro bably the fact that, as explained in IBerghofer et al] 
(1996), these authors computed the X-ray fluxes by using 
the count rates together with an energy conversion factor 
adopted on the basis of the hardness ratio. Though this 
method is certainly the best that can be done when the 
number of counts is low, it is inevitably less accurate 
than the present approac h . From Fig. |S| we however note 
that the IBerghofer et alJ l)l997f) relation obtained in the 
0.1-2.0 keV band for O-type stars remains qualitatively 
valid in the 0.5-10.0 keV band. Given the previous consid- 
erations, this agreement is probably a stroke of good fortune. 
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Figure 13. Same as Fig. |7| for O-type stars. Th e asterisks 
give in addition the lAlbacete Colombo et all J2003t) data (ob- 
tained in the 0.3-12.0 keV band) and the dashed-dotted line 
indicates our corresponding best-fit scaling law to these data 
(Eq. 1121 . The open star sym bols show two points that, follow- 
ing |Aftacete_Colombo_et_ay, were not included in the fit. 



earlier than 06, nor to bolometric luminosities larger than 
log Lboi = 39.5. Because of the limited range in Lboi (less 
than 1.5 order of magnitude), extrapolation towards earlier 
types should be considered with caution. Within these 
limitations and rejecting the peculiar sources (the two 
CWBs in the current sample), we note that no systematic 
discrepancy is observed for stars belonging to different 
luminosity classes. All the objects (except perhaps the late 
O-type dwarfs) seem to tightly follow the same canonical 
relation. This definitely suggests the existence of a common 
physical mechanism to generate the X-ray emission for 
'normal' (i.e. non peculiar) O-type stars. 



6.1.1 The natural dispersion around the canonical 
Z/x — Lboi relation 

From the ROSAT Al l -Sky survey data (IBerghofer et alJ 
1996), Bcrghof er et alJ |l$97?) derived a natural dispersion 
around the canonical relation of about 0.40 in the log — log 
plane (thus of about a factor 2.5 on the X-ray luminosity), 
a dispersion that was already much reduced compared to 
ISciortino et al] (ll990T) work. While part of this dispersion 
is probably due to the non-discrimination of the peculiar 
objects, it could also reflect the instrumental limitations 
of the time. Beyond the intrinsic quality of the data, both 
the homogeneity of the present sample and the larger value 
adopted for the low energy cut-off have probably played an 
important role in preserving a low level of dispersion and, in- 
deed, the dispersion observed around the canonical relation 
in NGC 6231 is definitely very limited. 

The current sample of O stars contains both single stars 
and binaries, and includes objects from different luminosity 
classes. However, it does not extend towards spectral types 

5 These values were computed using the HEASARC's 
W3PIMMS tool available at http://heasarc.gsfc.nasa.gov/ 
Tools/w3pimms .html. We used a Raymond-Smith model with 
a temperature of 1.0 keV, but the latter parameter has only a 
limited influence on the quoted results. 



Except for some possible observational artifacts (e.g. 
contamination by neighbouring sources, see Sect.0, the sole 
cause for the observed deviation from the canonical rela- 
tion in the present sample is X-ray emission from a wind 
interaction zone. This extra-emission however does not sig- 
nificantly affect all the O-type binaries. It is also the sole 
identified mechanism that, at our detection limit, yields a 
modulation of the X-ray flux in our star sample. In conse- 
quence, the present analysis suggests that the intrinsic X-ray 
emission from non-peculiar O-type stars does not show any 
significant variability and can be considered as constant for 
a given star. 

6.1.2 Constraints from the Carina region 

lAlbacete Colombo et all (I2003T) have recently presented a 
44 ks XMM-iVetufon observation of the Carina nebula field. 
As part of their analysis, they reported the canonical re- 
lation to be logLx = 1.07(±0.04) log L bo i - 6.2(±0.1) in 
the 0.3-12.0 keV domain. Unfortunately, the quoted rela- 
tion overestimates the X-ray luminosities compared to their 
actual measurements 6 . Using the data from their Table 6 
and considering the same objects as they did (11 out of 

6 Actually, it is possible that the quoted constant in the power- 
law is erroneous. 
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the 13 O stars located in the FOV), we have rederived the 
best-fit power-law relation. Applying an unweighted linear 
regression, one obtains : 

log L x = 1.02(±0.15) log L bo i - 6.9(±5.8) . (11) 

The typical dispersion around the relation is about 0.18 
in the log — log plane, thus very similar to the dispersion 
obtained from our data set. Adopting the scaling relation 
model, an unweighted fit yields 

log (Lx/L ho i) = -6.23 ±0.17 . (12) 

Comparing the latter residuals with those obtained from 
Eg. 1111 we obtained F x — 0.01. This clearly indicates that 
the power law does not provide any significant improve- 
ment of the fit compared to the scaling law. The latter 
relation is plotted in F ig. II 3 | together with the data of 
lAlbacete Colombo et all l)2003h . 

This new analysis of the lAlbacete Colombo et alJ (120031) 
data puts into light two important points. First, as in our 
case, the preferred log Lx — log Lboi relation takes the form of 
a sca ling law. This clearly contrasts with the previous s tudies 
(e.g. IChlehowski fc OarmamJliflflil iBerghofer et all Il997h 
which probably adopted the power law to provide an extra 
degree of freedom to the fit. The second point relies on the 
different nature of the considered samples. Our sample con- 
tains Q 6-Q9 stars belonging to d ifferent luminosity classes 
whereas lAlbacete Colombo et alJ <l2003t) data set is almost 
uniformly formed by main-sequence stars ranging from 03 
to 08.5. However, in both cases, the scaling law is preferred 
and both samples displayed a limited dispersion around the 
best-fit relation. We suggest that the slightly la rger disper- 
sion observed in lAlbacete Colombo et all i2003h data could 
come from a less uniform analysis and, in particular, from 
the handling of the correction for the ISM absorption. Ac- 
counting for the difference in the analysis (different spectral 
models used) and in the considered energy ranges, the two 
relations can be regarded as being in mutual agreement if 
about one fifth of the intrinsic X-ray flux is actually emitted 
in the 0.3-0.5 keV range. 

6.1. 3 On the origin of the Lx — Lboi scaling 

As pointed out by different authors already, the scaling be- 
tween Lx and Lboi could be an indirect effect. Indeed, in 
the framework of the embedded wind shock model, it is 
expected that the X-ray luminosity rather scales with the 
wind parameters or with some physica lly meaningful com- 
bination of them. ISciortino et alJ l|1990]) observed a correla- 
tion with the wind momentum (Muoo) and the wind kinetic 
luminosity (0.5Mw^,) but did not find a direct correlation 
with the sole mass-loss rate (M) nor with the terminal ve- 
locity (t>oo). It is only through the scaling of the properties 
of the radiatively driven winds with the stellar luminosity 
that t he dependence between Lx and L bo i is supposed to 
occur. iKudritzki et"afl dl99d) and lOwocki fc Cohenl Jl99fih 
later suggested that an important parameter could be the 
mean wind density, whi ch is related to M I twFrom the exo- 
spheric approximation, lOwocki fc Cohenl jl999l) established 
that the X-ray luminosity is expected to scale with the wind 
density parameter M/v x in the form Lx ~ {M/voo) 2 for 
optically thin winds and Lx ~ (M/«oo) for optically 



thick winds, with s being the index of the radial power law 
dependence of the X-ray fil li ng factor / ~ r s . 

Using the IVink et alJ fcOOl 1200 J) recipes, we esti- 
mated the density parameter (M/vaa), the wind momen- 
tum {Mvoa) and the wind kinetic luminosity (0.5Mu^) for 
the different O-type objects in our sample. We assumed 
i>oo = 2.6 w csc . Masses and temperatures we re taken from 
the ne w calibration for O star parameters of lMartins et alJ 
(2005) while the Lboi were taken from Sect.|H] Fig. lMl shows 
the resulting diagrams. Using the same sample of stars as 
for Eq. [5] we obtained: 

log L x = (0.402 ± 0.064) log (M/v x ) + (27.46 ± 0.71) (13) 
log L x = (0.513 ± 0.075) log (Mv x ) + (17.74 ± 2.07) (14) 

log L x = (0.574 ± 0.084) log (o.5Mi&) + (11.45 ± 3.01) (15) 

with a typical residual dispersion of 0.16. Eg. 1131 indicates a 
filling factor index s ~ —0.6. 

Wind te rminal velocities and sta r mass-loss rates were 
measured bv lHowarth fc Prinial il989T) for five of the bright- 
est stars of our sample. The corresponding observational 
wind parameters are displayed in Fig. 1141 The agreement 
with the derived relation is rather acceptable although slight 
systematic upward shifts are observed. Only HD 152248 
shows strong deviations from these relations, but the deriva- 
tions of M and Voo have probably been biased by the ongoing 
wind-wind interaction. This clearly supports the idea that 
the O star X-ray luminosity also scales with various combi- 
nations of the wind parameters M and Voc . 

However, we emphasize that the canonical relation ex- 
pressed in Eq. |2]and the Eos. H3l to ll5l are not independent 
but rather reflect a different point of view on the scaling 
of the X-ray luminosities with the O-type star properties. 
Indeed the mass-loss recipes of IVink et alJ analytically ex- 
presses M in terms of Lboi, T e g and the stellar mass M. 
Using Eq. [5] to express Lboi in terms of Lx in the mass 
loss recipes, we re-derived analytically very similar slopes 
as those quoted in Eqs. 1131 to 1151 For example, subtracting 
log Woo from logM almost removes the first order depen- 
dency in T c ff. The remaining relation is dominated by the 
Lboi term. From these considerations, it appears that one 
can naturally convert the canonical relation between the X- 
ray and bolometric luminosities into a relation linking the 
X-ray luminosity to the wind paramet ers. For example, Lx 
may scales with M/v-a as proposed bv lOwocki fc Cohenl on 
the basis of more physical considerations. 

It is beyond the scope of this paper to pursue the quest 
for the physical origin of the scaling between the X-ray and 
the bolometric luminosities. Though the physical processes 
responsible of the X-ray emission in hot stars are not directly 
related to Lboi, the scaling between these two parameters is 
an observational fact. Compared to parameters such as M 
and Voa, or any combination of these, the bolometric lumi- 
nosity is probably much easier to measure and can be deter- 
mined with a better accuracy. We note that this relation is 
now tightly constrained and still awaits to be explained on 
firm theoretical grounds. 
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Figure 14. log Lx plotted ve rsus differen t wind param eters. The filled symbols have the same meanings as in Fig.||J Their abscissa were 
computed using the recipes of IVink et alj (2000, 200l| ), Dashed lines show the obtained best- fit relations (Eqs. IT51 to 1151 . The asterisks 
were computed using the M and v x measurements of lHowarth &; Prinial il9S9f) for the available objects. 



6.2 The B-type X-ray emitters 

Among the ~90 B-type stars in the XMM-Newton FOV of 
the present campaign, less than 20% could be associated 
with an X-ray counterpart. However, except for the earliest 
sub-types, the B stars are not supposed to emit in the X-ray 
domain. Their stellar winds are much weaker than those of 
O-type stars and should thus not contribute to the X-ray 
emission. In addition, no coronal emission is expected as 
these stars are not supposed to have the required convective 
zones. From our data, either the B-type stars are intrinsic 
X-ray emitters or the detected X-ray emission is actually 
associated with a later-type or PMS object along the same 
line of sight or physically linked to the B star. The detected 
dependence between the B star bolometric luminosities and 
the associated X-ray luminosities, with a linear-correlation 
coefficient of ~0.75 (thus significant at the 0.01 level) is in 
favour of the first hypothesis. On the other hand, the fact 
that we only detect a fraction of the B stars at a given spec- 
tral type is clearly against. Indeed, if the X-ray luminosity 
was actually an intrinsic property of the star, it should be 
expected that similar stars would display similar Lx- 

In regard of these two hypotheses, the study of the vari- 
ability of the individual objects presented in Sect. |S] is par- 
ticularly relevant. It indicates that, among the 11 B-type 
stars used to derived the Lx — Lboi relations, four are dis- 
playing a flaring activity while two others present long term 
variations. While flares h ave been recently obser ved on a 
magnetic B-type star fe.g. lGroote fc Schmitdl2004f) . this re- 
mains rather exceptional and such a behaviour is more likely 
related to the presence of a PMS companion, either in a bi- 
nary system or located along the line of sight. Indeed, the 
X-ray spectra of the "B-type" sources in our sample are 
strikingly similar to those of the numerous PMS stars in the 
FOV (see our forthcoming analysis for more details about 
the PMS sources). 

As shown in lPaperll only about 40 X-ray emitters over 
the 610 sources are expected to be foreground objects. The 
possibility for such an object to be located by chance on the 
same line of sight of a B-type star is thus very limited and 
can not explain the number of associations between an X-ray 
source and a B-type star. The density of PMS stars within 



the Sco OB 1/NGC 6231 complex is much larger. However, 
the limited cross-correlation radius adopted ensures a very 
low probability (< 0.03) that such an X-ray emitting PMS 
star be located by chance on the same line of sight. There- 
fore, the hidden binary companion becomes a privileged ex- 
planation. 

While about 50% of the B-type stars are expected to ac- 
tually be in binary or multiple systems, typical X-ray emis- 
sion from a late-type companion is probably too faint to be 
detected at the distance of NGC6231. Hence, only the sys- 
tems hosting a PMS companion are expected to be seen in 
the X-ray domain. This could explain the difference between 
the expected fraction of binaries and the observed detection 
rate among B stars. In this regar d , our conclusions some- 
what rejoin those of IStelzer et al] i2005t) . who stated that 
the X-ray emission of weak wind stars is mostly consistent 
with magnetic activity from known or unseen low-mass com- 
panions. 

Though the above discussion is clearly in favour of the 
PMS companion scenario, the observed correlation, at a sig- 
nificance level of 0.01, between the bolometric luminosity of 
the B star and the associated X-ray luminosity is intriguing. 
A possible scenario could imply some connection between 
the B-type star and its PMS companion. Though highly pu- 
tative, such a mechanism could justify the observed link be- 
tween the B-type properties (Lboi) and the observed X-ray 
emission from the PMS companion. However, the observed 
relation might also result from a border effect. Indeed, if the 
real distribution of the emission 'associated' with B stars 
has about the form of a cluster of points and if one only 
observes the upper part of the distribution, one might po- 
tentially detect an apparent relation between the two con- 
sidered variables although, actually, it results from an ob- 
servational cut-off. 



6.3 The demarcation line 

In this paper, we have voluntarily restrained our analysis to 
objects whose spectra were of sufficient quality to be fitted 
with 2-T models. This actually led to the rejection of the 
early B-type stars that were only marginally detected. It 
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is not impossible that some of these actually form the low 
luminosity tail of the Lx — Lboi relationship derived for O- 
type stars. Another possibility is that, in this tail, the stellar 
winds are becoming so weak that only weak shocks are pro- 
duced, yielding thus a possibly softer emission (see e.g. the 
case of HD 152235 in Sect. [^J . Our approach was clearly not 
designed to investigate the location of the demarcation line. 
Still, the separation between the two different behaviours 
clearly occurs at Lboi ~ 10 38 ergs -1 . Furthermore, this value 
clearly separates the O and B-type stars of our sample. 



7 CONCLUSIONS 

In this second paper of the series, we have pursued the 
analysis of the X-ray data con cerning the young open clus- 
ter NGC6231. While IPaper il focused mainly on the detec- 
tion and identification of the numerous X-ray sources in the 
XyiM-Newton FOV, this paper was devoted to the proper- 
ties of the rich early-type star population. A detailed cen- 
sus of the OB-type stars within the FOV resulted in more 
than one hundred objects identified. Using a limited cross- 
correlation radius of 2'.'5, about one third of them could be 
associated with an X-ray counterpart. Among these, the 15 
O-type stars/binaries are all detected in the X-rays as soft 
and usually bright sources characterized by mekal temper- 
atures of kT = 0.3 and 0.7 keV. On the other hand, the 
B-type star detection rate only amounts to about 20%. Com- 
pared to the O-type stars, the B-type stars have a similar 
low energy component but their second temperature is well 
above 1 keV. The B-type X-ray emitters are thus signifi- 
cantly harder than the O-type sources. 

O- and B-type emitters clearly present different be- 
haviours in the log Lx — log Lboi diagram, though both types 
draw up a linear relation in the log-log plane. The separation 
between the two sub-sets is located at about log Lboi = 3 8 
(ergs -1 ), as previously suggested bv lBerghofer et al.Nl997l) . 
The dispersion around the expected linear relation is ap- 
parently quite small. In the O-type star sample, the two 
objects that show the largest deviations are known to dis- 
play an extra -emission compon ent due to a wind interaction 
(HP 152248. ISana et all 12004 CPD -41° 7742. ISana etafl 
2005). These were thus excluded from the subsequent dis- 
cussion. We showed that, for the O-type stars, the X-ray 
luminosities are scaling with the bolometric luminosities. In 
the 0.5-10.0 keV energy range, we obtained: 

log Lx - log Lboi = -6.912 ± 0.153. 

We also found that a power law relation did not provide 
any significant improvement to the quality of the fit. The 
obtained dispersion around this new canonical relation is 
very limited. It becomes even smaller when excluding the 
'cooler' (i.e. with a spectral type later than 09) O dwarfs 
from the fit. In this case, the typical dispersion drops to 
about 0.087 in the log — log plane, thus corresponding to 
only 20% on the X-ray luminosities. Within our sample, the 
only identified mechanism that provides a significant devia- 
tion from this relation is extra-emission produced in a wind 
interaction region. It is also the sole mechanism that, at 
our detection threshold, produces a significant variability of 
the observed fluxes in our O-type star sample (see Sect.0. 
Though relatively limited, the present sample suggests thus 



that the intrinsic X-ray emission from O-type stars is very 
tightly correlated with their bolometric luminosity. Beyond 
the two strong CWB systems, our sample is formed by single 
stars and binaries belonging to different luminosity classes. 
Though not extending towards sub-types earlier than 06, 
they all seem to follow the canonical relation, suggesting 
thus a common mechanism for X-ray production in these 
objects. 

We also provide a new analysis of recent 
flux measurements ob t ained in the Carina region 
iAlbacete Colombo et all l2003t) . This sample is mainly 
formed by main-sequence stars ranging from 03 to 08.5. 
We note that this new analysis confirms much of our 
present conclusions. The best Lx — Lboi relation is indeed 
in the form of a scaling law rather than a power law. The 
dispersion around the obtained relation is very limited 
and the difference with our own relation might possibly be 
accounted for by the different energy ranges considered. 

We emphasize that this apparent scaling might indi- 
rectly result from a scaling of the X-ray luminosity with 
the wind properties, themselves scaling with the bolometric 
luminosities for these stars with radiati vely driven winds . 
Clues for this are provided by the work of lVink et alJ (l2000l . 
2001) which has allowed us to naturally convert the scal- 
ing of the X-ray luminosity with the bolometric luminosity 
into a scaling with the wind parameters. Nonetheless the 
Lx — Lboi relation probably remains the most accurate ob- 
servational constraint to link the intrinsic X-ray emission of 
the O-type stars with their fundamental properties. 

Turning to B-type stars, the fact that only about one 
quarter of the stars of a given spectral sub-type are actually 
associated with an X-ray source argues strongly against X- 
ray emission being an intrinsic property of these stars. We 
however note that we still observed a linear relation that 
links logLx and log Lboi: 

log L x = (0.22 ± 0.06) log L bo i + 22.8(±2.4) 

The dispersion around this relation is quite limited (~0.14) 
and the linear pattern is also seen in the different energy 
sub-ranges considered. From our analysis, the most probable 
explanation points towards the X-ray emission originating 
from a low mass PMS physical companion. The observed 
relation between Lx and Lboi remains however a puzzle 
and could eventually result from an observational effect. 
Alternatively it could be linked to some particularities of 
the B-type stars in NGC 6231 or to a putative interac- 
tion between the PMS object and its B-type companion, 
yielding an X-ray emission partly governed by the intrinsic 
properties of the B-type primary. 

Finally, we note that the separation line between the 
O- and "B-type" behaviours (around Lboi = 10 38 ergs -1 ) is 
poorly mapped by our present sample. The adopted criterion 
to preserve the homogeneity of the spectral analysis (requir- 
ing at least 2-T models) has led to the rejection of the objects 
in the transition zone. The extent of the canonical relation 
towards lower luminosities probably deserves a more partic- 
ular attention. Dedicated observations, combined with the 
already observed fields, could help to increase the number of 
objects in this zone. This is probably a necessary condition 
to answer with more details this still open question. 



20 H. Sana et al. 



ACKNOWLEDGMENTS 

The authors are greatly indebted to the 'Fonds National 
de la Recherche Scientifique' (FNRS), Belgium, for multi- 
ple supports. This research is supported in part by contract 
P5/36 "Pole d'Attraction Interuniversitaire" (Belgian Fed- 
eral Science Policy Office) and through the PRODEX XMM 
and INTEGRAL contracts. The SIMBAD and WEBDA 
databases and the Vizier catalogue access tool (CDS, Stras- 
bourg, France) have been consulted for the bibliography and 
for the purpose of object cross-identification. 



REFERENCES 

Albacete Colombo J. F., Mendez M., Morrell N. I., 2003, 
MNRAS, 346, 704 

Arentoft T., Sterken C, Knudsen M. R., Freyhammer 
L. M., Duerbeck H. W., Pompei E., Delahodde C. E., 
Clasen J. W., 2001, A&A, 380, 599 

Arnaud K. A., 1996, in Jacoby G, Barnes J., eds, Astro- 
nomical Data Analysis Software and Systems V Vol. 101 
of ASP Conf. Ser., XSPEC: The First Ten Years, p. 17 

Balona L. A., 1983, MNRAS, 203, 1041 

Balona L. A., Engelbrecht C. A., 1985, MNRAS, 212, 889 

Balona L. A., Shobbrook R. R., 1983, MNRAS, 205, 309 

Baume G, Vazquez R. A., Feinstein A., 1999, A&AS, 137, 
233 

Berghofer T. W., Schmitt J. H. M. M., Cassinelli J. P., 

1996, A&AS, 118, 481 (Erratum: 1997, A&AS, 121, 212) 
Berghofer T. W., Schmitt J. H. M. M., Danner R., 

Cassinelli J. P., 1997, A&A, 322, 167 
Bevington P., 1969, Data Reduction and Error Analysis for 

the Physical Sciences. McGraw-Hill Book Company, Inc., 

USA 

Bohlin R. C, Savage B. D, Drake J. F., 1978, ApJ, 224, 
132 

Braes L. L. E., 1967, Bulletin of the Astronomical Institute 

of the Netherlands Supplement Series, 2, 1 
Cassinelli J. P., Olson G. L., 1979, ApJ, 229, 304 
Cassinelli J. P., Waldron W. L., Sanders W. T., Harnden 

F. R., Rosner R., Vaiana G. S., 1981, ApJ, 250, 677 
Cassinelli J. P., Cohen D. H., MacFarlane J. J., Sanders 

W. T., Welsh B. Y., 1994, ApJ, 421, 705 
Chlebowski T., Garmany C. D, 1991, ApJ, 368, 241 
Chlebowski T., Harnden F. R., Sciortino S., 1989, ApJ, 

341, 427 

Cohen D. H., Cassinelli J. P., MacFarlane J. J., 1997, ApJ, 
487, 867 

Cohen D. H., de Messieres G. E., MacFarlane J. J., Miller 
N. A., Cassinelli J. P., Owocki S. P., Liedahl D. A., 2003, 
ApJ, 586, 495 

Feigelson E. D, Montmerle T., 1999, ARA&A, 37, 363 
Feldmeier A., Puis J., Pauldrach A. W. A., 1997, A&A, 
322, 878 

Gagne M., Oksala M. E., Cohen D. H., Tonnesen S. K., ud- 
Doula A., Owocki S. P., Townsend R. H. D, MacFarlane 
J. J., 2005, ApJ, 628, 986 
Garcia B., Mermilliod J. C, 2001, A&A, 368, 122 
Garrison R. F., Schild R. E., 1979, AJ, 84, 1020 
Groote D., Schmitt J. H. M. M., 2004, A&A, 418, 235 



Harnden F. R., Branduardi G, Gorenstein P., Grindlay J., 
Rosner R., Topka K., Elvis M., Pye J. P., Vaiana G. S., 
1979, ApJ, 234, L51 
Hill G, Crawford D. L., Barnes J. V., 1974, AJ, 79, 1271 
Houk N., 1978, Michigan catalogue of two-dimensional 
spectral types for the HD stars. Ann Arbor : Dept. of 
Astronomy, University of Michigan : distributed by Uni- 
versity Microfilms International, 1978 
Howarth I. D., Prinja R. K., 1989, ApJS, 69, 527 
Howk J. C, Cassinelli J. P., Bjorkman J. E., Lamers 

H. J. G. L. M., 2000, ApJ, 534, 348 
Isobe T., Feigelson E. D., Nelson P. I., 1986, ApJ, 306, 490 
Kaastra J., 1992, An X-Ray Spectral Code for Optically 
Thin Plasmas. (Internal SRON-Leiden Report, updated 
version 2.0) 

Kramer R. H., Cohen D. H., Owocki S. P., 2003, ApJ, 592, 
532 

Kudritzki R. P., Palsa R., Feldmeier A., Puis J., Pauldrach 
A. W. A., 1996, in Roentgenstrahlung from the Universe, 
Eds. H.U. Zimmermann et al., MPE Report 263, p. 9 
Laval A., 1972, A&A, 21, 271 
Levato H., Malaroda S., 1980, PASP, 92, 323 
Levato H., Morrell N., 1983, Astrophys. Lett., 23, 183 
Long K. S., White R. L., 1980, ApJ, 239, L65 
Lucy L. B., 1982, ApJ, 255, 286 

Martins F., Schaerer D., Hillier D. J., 2005, A&A, 436, 1049 
Mayer P., Harmanec P., Lorenz R., Drechsel H., Eenens P., 

Corral L. J., Morrell N., 2001, in The Influence of Binaries 

on Stellar Population Studies, ASSL Vol. 264, p. 567 
Mewe R., Gronenschild E. H. B. M., van den Oord G. H. J., 

1985, A&AS, 62, 197 
Owocki S. P., Cohen D. H., 1999, ApJ, 520, 833 
Owocki S. P., Castor J. I., Rybicki G. B., 1988, ApJ, 335, 

914 

Pallavicini R., Golub L., Rosner R., Vaiana G. S., Ayres 

T., Linsky J. L., 1981, ApJ, 248, 279 
Perry C. L., Hill G, Younger P. F., Barnes J. V., 1990, 

A&AS, 86, 415 
Perry C. L., Hill G, Christodoulou D. M., 1991, A&AS, 

90, 195 

Porquet D., Mewe R., Dubau J., Raassen A. J. J., Kaastra 

J. S., 2001, A&A, 376, 1113 
Raboud D., 1996, A&A, 315, 384 
Reed B. C, 2003, AJ, 125, 2531 

Rosner R., Golub L., Vaiana G. S., 1985, ARA&A, 23, 413 
Rufener F., Bartholdi P., 1982, A&AS, 48, 503 
Sana H., 2005, PhD thesis, Liege University 
Sana H., Rauw G, Gosset E., 2001, A&A, 370, 121 
Sana H., Hensberge H., Rauw G, Gosset E., 2003, A&A, 
405, 1063 

Sana H., Stevens I. R., Gosset E., Rauw G, Vreux J.-M., 

2004, MNRAS, 350, 809 
Sana H., Antokhina E., Royer P., Manfroid J., Gosset E., 

Rauw G, Vreux J.-M., 2005, A&A, 441, 213 
Sana H., Naze Y., Gosset E., Rauw G, Sung H., Vreux J.- 

M., 2006a, in Moffat A., St-Louis N., eds, Massive Stars 

in Interacting Binaries ASP Conf. Ser., pp 5p., in press 
Sana H., Gosset E., Rauw G, 2006b, MNRAS, in press 
Sana H., Gosset E., Rauw G, Sung H., Vreux J.-M., 2006c, 

A&A, 454, 1047 (Paper I) 
Schild R. E., Neugebauer G, Westphal J. A., 1971, AJ, 76, 

237 



An XMM-Newton view of NGC6231. 



Schmidt-Kaler T., 1982, Physical Parameters of the Stars. 
Vol. 2b of Landolt-Bornstein, Numerical Data and Func- 
tional Relationships in Science and Technology, New Se- 
ries, Group VI, Springer- Verlag, Berlin 

Schmitt J. H. M. M., Golub L., Harnden F. R., Maxson 
C. W., Rosner R., Vaiana G. S., 1985, ApJ, 290, 307 

Schulz N. S., Canizares C, Huenemoerder D., Tibbets K., 
2003, ApJ, 595, 365 

Sciortino S., Vaiana G. S., Harnden F. R., Ramella M., 
Morossi C, Rosner R., Schmitt J. H. M. M., 1990, ApJ, 

361, 621 

Seward F. D., Chlebowski T., 1982, ApJ, 256, 530 
Seward F. D., Forman W. R., Giacconi R., Griffiths R. E., 

Harnden F. R., Jones C, Pye J. P., 1979, ApJ, 234, L55 
Stelzer B., Huelamo N., Hubrig S., Zinnecker H., Micela 

G., 2003, A&A, 407, 1067 
Stelzer B., Flaccomio E., Montmerle T., Micela G., 

Sciortino S., Favata F., Preibisch T., Feigelson E. D., 2005, 

ApJS, 160, 557 
Stickland D. J., Lloyd C, Penny L. R., 1997, The Obser- 
vatory, 117, 213 
Sung H., Bessell M. S., Lee S., 1998, AJ, 115, 734 
ud-Doula A., Owocki S. P., 2002, ApJ, 576, 413 
Vaiana G. S., Cassinelli J. P., Fabbiano G., Giacconi R., 

Golub L., Gorenstein P., Haisch B. M., Harnden F. R., et 

al, 1981, ApJ, 245, 163 
Vink J. S., de Koter A., Lamers H. J. G. L. M., 2000, A&A, 

362, 295 

Vink J. S., de Koter A., Lamers H. J. G. L. M., 2001, A&A, 
369, 574 



This paper has been typeset from a TfjjX/ LTgX file prepared 
by the author. 



